UNCLASSIFIED 


AD  NUMBER 


AD251478 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release^  distribution 
unlimited 


EROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  SEP  1960. 
Other  requests  shall  be  referred  to  Air 
Force  Special  Weapons  Center^  Kirtland 
AFB,  NM. 


AUTHORITY 


AFWL  Itr,  25  May  1970 


THIS  PAGE  IS  UNCLASSIEIED 


UNCLASSIFIED 


in  251  478 


fieo'ioduced 
ixf  ike 


ARMED  SERVICES  TECHNICAL  INFCRMAHON  AGENCY 
ARLINGTON  HALL  STAT1(M« 

ARLINCra  12,  VIRGINU 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  In  connection  with  a  definitely  related 
govertu5'‘nt  procurement  operation,  the  U.  S. 
Govemiaent  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  In  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  Is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  bolder  or  any 
other  person  or  corporation,  or  conveying  any  ri^ts 
or  pemlssion  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  einy  way  be  related 
thereto. 


AFSWC-TN-60-36 


SWC 

TN 

60-36 


HEADQUARTERS 


//  ■ 


111  ItSEliei  UO  DIVELOPaEIT  COHMIl 

u,.  ' 

KIITLIINI  Ml  FORCE  USE,  NEW  MEXICO 


^  Technical  Note  ,  , 

DESIGN  AND  ANALYSIS  OF  FOUNDATIONS 
'  FOR  PROTEOT’*  St-lTuCTUAES  ■  K 

Inter' in  Tei  Imical  Report 
by 

K.  E.  McKee 

Armour  Research  Foundation  of 
Illinois  Institute  of  Technology 
Technology  Co»iter 
Chicago  16,  I’*  ois 


September  1960 


AFSVfC-TN-60-36 


DESIGN  AND  ANALYSIS  OF  FOUNDATIONS 
FOR  PFCOTECTIVE  STRUCTURES 
Interim  Technical  ReporJ: 

by 

K.  E.  McKee 

Armour  Research  Foundation  of 
Illinois  Institute  of  Technology 
Technology  Center 
Chicago  16,  Illinois 

September  196  Q 


Research  Directorate 
AIR  FORCE  SPECIAL  WEAPONS  CENTER 
Air  Research  and  Development  Command 
Kirtland  Air  Force  Base,  New  Mexico 


Approved; 


Project  No,  lOS’O 

Task  No.  10803 

Contract  No.  AF  29(601)-2561 

ARF  Project  No.  K193 


PER^  L.  HUIE 
Col^l  USAF 
Chief,  Structures  Division 


LklONARDA.  EDD’i 
Colonel  USAF 

Director,  Research  Directorate 


TN-60-36 


ABSTRACT 

The  behavior  of  footings  subjected  to  time- dependent  forces  has  been 
the  subject  of  continuing  research.  The  ultimate  bearing  capacity  under  such 
loading  conditions  and  the  dynamic  behavior  beyond  this  ultimate  capacity 
arf  both  of  interest.  An  attempt  has  been  made  in  the  subject  Investigation 
to  combine  and  correlate  laboratory  experiments  with  theoretical  studies. 

Two-  and  three-dimensional  experiments  have  been  conducted  on 
small  footings  in  the  laboratory  to  observe  their  behavior  and  to  obtain 
quantitative  information.  An  apparatus  was  developed  for  applying  dynamic 
forces  to  the  footings.  This  apparatus,  which  is  relatively  simple,  has  made 
possible  the  application  of  loads  having  various  rise  times,  decays,  and 
durations.  Force-time  and  displacement-time  records  have  been  obtained 
in  forms  suitable  for  analysis,  and  Fastax  movies  of  footings  failing  under 
dynamic  loads  have  also  been  taken. 

The  behavior  of  footings  subjected  to  dynamic  loads  has  been  studied 
analytically  also.  The  possibility  of  applying  the  plasticity  theory  or 
limit  anaiyoia  has  been  considered.  Other  loadings  and  various  failure  modes 
also  have  been  investigated.  Use  %\dll  be  made  of  the  experimental  data  in 
conjunction  with  this  theoretical  work  during  the  remainder  of  the  program. 
Based  on  this  work,  additional  experimental  and/or  theoretical  research  will 
be  conducted  as  required. 


PUBLICATION  REVIE'Y 
This  report  has  been  reviewed  and  is  approved. 


PHEFACE 


This  is  the  interlJi  technioal  report  on  Contract  ITo.  AF29(&)1)-2561,  . 
Project  i080|  "Design  and  Analysis  of  Foundations  for  Protective  Strueiairas", 

The  object:  j  of  this  research  progran  is  to  Investigate  the  problens  associated 
with  the  design  and  analysis  of  foundations  for  protective  structures  which 
are  subjected  to  dynanio  loads  from  nuclear  blast.  The  current  project,  which 
was  initiated  in  February,  196C,  to  a  large  extent  is  a  continuation  of  ressaroh 
done  on  an  earlier  oontraot,  AF29(601)-ll6l,  which  h^d  the  sane  title. 

This  interim  technical  report  covers  the  work  done  on  Ihe  project 
during  the  seven  months  following  the  inception  of  the  con-tract.  In  organizing 
this  report,  an  at-tempt  was  made  to  collect  the  pri.mary  technical  results  into 
appendices  which  subsequently  can  be  incorporated  into  the  final  report. 

Personnel  who  have  contributed  to  the  work  covered  in  this  report 
Includai  R,  L,  Ghiapetta,  C.  J.  Costantlno,  P,  0,  Hodge,  A,  Humphreys,  K,  E, ,  McKee 
R,  D,  Rowe,  R.  W,  Sauer,  S.  T.  Selig,  S,  Shenknan,  and  E,  Vey.  Credit  should 
also  be  given  to  Mr,  C,  Vlehle  and  Mr,  H,  Mason  of  AFSWC  for  their  criticisms 
and  suggestions  which  have  aided  this  project, 


ARMOUR  RESEARCH  POUNDATIQN  JLLINOIS  IMSTITUTE  OP  TECHSOLOGY 


iv 


ARF  Project  No,  8193-7 
Interim  Report 


TABLE  OF  CONTENTS 


ABSTRACT  . . . . . 

PREFACE  . . . . 

Chapter 

1  UrlRODUCTION . 

A.  Qenerol  . . 

B.  Ob . . . . . 

C.  Teshnlaal  Approach . . 

D.  Dieouseion  of  the  Problem . . 

El  Report  Crganlzatlon* •  I . . 

2  r.-lBORETICAL  Ur/ESTIOATION . 

A.  General • ....  n  i  n  .ki  .  m n  . 

B.  Plaatioity  Theory  in  Soil  Meohanios 

C.  Dynamic  Response I . 

D.  Non«Vartioal  Loads . . . . 

3  SXPEPJMEJrrAL  UnrSSTIQATrON . . 

A.  General . . . . 

B.  Material  Propei^ies . * . 

C.  Tvo-Dimensional  Studies . . . 

D.  Three-Diinensional  Studies . . 

li  CONCLUSIONS . . 


iv 


1 

«L 

1 

2 

2 

3 

u 

li 
h 
■  5 
8 

12 

12 

12 

13 

13 

Ih 


ARMOUR  RESEARCH  C O U N  D  A  T  I  O N  DC  ILLINOIS  :  N  S  T  I  T U r  E  OE  TECHNOLOGY 


V 


ARE  Proj-eot  No.  8193-7 
Interim  Report 


TABLE  OF  CONTENTS,  cont. 

Appendix 

A  APPLICATION  OF  PLASTICITY  AND  LIMIT  ANALYSIS  TO 

FOOTINO  DESICaJ  . . . . . . . .  16 

Atl  Introduction. . . . . .  17 

A«2  Plasticity  Theory..... . .  16 

A>3  Limit  Theories...., . . . .  21 

A.ii  Previous  Studies....... . . . .  2$ 

k,$  Illustrative  Problem................. ............  28 

A. 6  Dynamic  Problem................ . . .  32 

A.  7  References... . . . . . . .  33 

B  THE  METHOD  OF  CHARACTERISTICS  APPLIED  TO  PRC3L2<S 

OF  PLANE  PLASTIC  STRAIN  ’rflTH  INStTIAL  EFFECTS., . .  35 

B. l  Introduction.... . . . 36'"" 

B.2  Basic  Equations,... . . .  37 

B.3  Differerwe  Equations  for  Small  Deformations*..,,,  ’  U2 

B.iiExasqsle............... ................ ...........  1^8 

B.5  Conclusions. . . »,,,■-$! 

B, 6,  References.  32 

C  FOOTDMS  SUBJECTED  TO  ECCENTRIC  LOADS . 38 

C. l '  Introduction, . . . 3^ 

^-^,2  ’  Static  Loads. ....................................  60 

''Dynatnic  Loads.. . . 6l 

D  PROV':;;iriES  OF  OTTAWA  SAND . .  66- ' 

E  Tl’/O-DIMENSTONAL  STUDIES . 73 

E.l  Introduction. . 7h 

E,2  Static  Tests...... . . . .  7h 

E. 3  lynamic  Tests . . . 73 

F  THREE-DIMENSIONAL  STUDIES .  88 

F. l  Introduction . . . 89 

F.2  Static  Tests....,......*. . . . .  89 

F,3  Dyi^ic  Tests . . . . .  90 

G  DYNAMIC  LOADING  APPARATUS . 92 


armour  research  poundation  op  Illinois  institute  op  technologt 

vi  ARF  Project  No.  8193-7 

Interim  Report 


LIST  OF  ILLUSTHATIONS 


Figure  Paga 


1 

SIMPLIFIED  IWO-SIDED  FAILURE  MODS  BASED  ON 

PRANDTL '  S  SOLUTION. . . . 

10 

2 

SIMPLIFIED  TkTO-SIDED  FAILURE  MODE  BASED  ON 

HILL'S  SOLUTION . 

11 

APPENDICES 

A.l 

AN  EQUILIBRIUM  SOLUTION . . . 

3h 

A.2 

A  VELOCITY  SOLUTION., . . . 

3k 

B,1 

CUHVILTNSAR  COORDINATE  SYSTSi . . 

55 

B.2 

TYPICAL  LIESH  IN  CHARACTERISTIC  NET . 

55 

B.3 

PUNCH  ON  HALF- SPACE <  DIMENSIONS  AND  STATIC 
SOLUTION . . . . 

56 

B,h 

PUNCH  ON  HALF-SPACE j  NATURE  OF  SjLUTION  AND 
BOUNDARY  CONDITIONS . . 

57 

C.l 

ONE-SIDED  FAILURE . . . 

63 

C.2 

STATIC  LOAD  vs.  FOOTING  WIDTH  AS  A  FUNCTION 

OF  8  . . 

6h 

C.3 

RADIUS  vs.  FOOTHra  WIDTH  AS  A  mCTION  OF  e 

,65' 

D.l 

GRAIN  SIZE  DISTRIBUTION  FOR  OTTAWA  SAIH),  BOTH 

NEW  SAND  AND  ORIGINAL  SAND . . . 

70 

D.2 

FAILURE  ENVELOPES  FOR  OTTAWA  SAND  AS  DETERMINED 
FROM  direct  SHEAR  TESTS .  • 

71 

D.3 

RESULTS  OF  TRIAXIAL  TF.STS  ON  OTTAWA  SAND. . . . 

72 

E.l 

FOOTING  TESTS  IN  GLASS  BOX,  DENSE  SAND,  3'-in. 

WIDE  Fooirao . 

77 

A  B  M  0  U  6  StSCABCH  BOUNOATlON  OC  ILLJNOIS  INSTITUTE  OF  T  C  C  Ft-N  O  L  O  G  T 


vil 


m  Project  No.  8193-7 
Interim  Report 


LIST  OF  ILLUSTRATIONS,  cont. 


E.2  FOOTING  TESTS  IN  GLASS  BOX,  LOOSE  SAND,  3-in. 

WIDE  FOOTING . . .  78 

E.3  STATIC- FOOTING  FAILURE  ON  DENSE  SAND . 79 

E.3  STATIC  FOOTING  FAILURE  ON  DENSE  SAND .  80 

E.ii  STATIC  FOOTING  FAILURE  ON  LOOSE  SAND . .  ,  8l 

E.I4  ..STAITO  FOOTING  FAILURE  ON  LOOSE  SAND.... .  82 

E,5  FOOTING  BmVIOR  ON  COHESIVE  SOIL .  83 

E.6  COHESIVE  SOIL  IN  GUSS  BOX  BEFORE  TEST . 81j 

E.?  COHESIVE  SOIL  IN  GLASS  BOX  DURING  STATIC  TEST,...  81i 

E.8  COHESIVE  SOIL  AFTER  GUSS  IN  BOX  FAILED .  8i* 

E.';  DYNAMO  TEST  OF  TOOTING  ON  DENJi;  OTTAwA  LL'JJD .  ' 

E.9  DYNAMIC  TESJ  OF  FOOTING  ON  oaSK  OTLl./A  Sa.'.D .  80 

E. IO  DYN*U>4lC  TEST  OF  FOOTING  ON  LOOO^  OTTA/iA  LLi.D .  87 

F. l  li-in.  SQUARE  FOOTING  “AFTER  FATLURi-; .  91 

F. 2  3-in.  SQUARE  FOOTING  AFTER  FAILURE .  91 

G, 1  SKETCH  OF  DYNAMIC  APPARATUS . 95 

0.2  DYNAMIC  APPARATUS  ON  GUSS  BOX. . 95 

G .3  DYNAMIC  APPARATUS  ON  SAND  BOX. . 96 

G.Ii  RECORDS  FOR  U-in.  x  Mn.  FOOTING,  TEST  NO.  6,...  97 

0,5  RECORDS  FOR  L-in.x  l4-tn,  FOOTING,  TEST  NO.  2 .  93 

0.6  TYPICAL  CONSOLIDATED  RLCOHDS,  TEST  NO,  U  99 

G.7  TYPICAL  CONSOLTDATED  RECORDS,  TEST  NO.  10 .  100 


A8M0UB  RESEARCH  POUNOATl'-ON  OP  rUtlNOIS  INSTITUTE  OP  TECHN 


O  L  O  G  T 


viii 


ARF  Project  No.  3193-7 
Interim  Report 


Chapter  1 
INTRODDCTION 


A.  General 

The  purpose  of  this  research  program  is  to  (insider  °the  design 
and  analysis  of  foundation  for  protective  structures  subjected  to  dynamic 
loads  from  nuclear  blast".  Witiiin  this  general  research  area,  thai*e  are  almost 
unlimited  technical  problems  that  could  be  considered.  This  report  is  concerned 
with  the  wor^  done  since  the  inception  of  the  present  program  in  February,  I960. 

In  this  report,  it  is  assumed  that  the  reader  is  familiar  vLth  the 
re^^lts  of  the  previous  research  conducted  by  ARF  in  the  same  problem  area 
(Contract  AF29(60lVll6l)  For  all  practical  purposes,  the  research  reported 
herein  represents  an  extension  of  the  previous  program.  Despite  tills  oloae 
association,  it  has  been  the  aim  to  make  this  report  readable  without  supple¬ 
mentary  references. 

This  report  is  organized  to  present  a  self-contained  progress 
report  of  the  technical  work.  To  the  maximum  extent  possible,  detailed 
technical  presentations  are  confined  to  the  appendices  with  the  boc^y  of  the 
report  limited  to  general  presentations,  conclusions,  eto.  It  is  to  be  expected 
that  the  reader  will  mako  use  of  thee©  appendices,  and  little  technical  infor¬ 
mation  is  reproduced  in  the  main  body. 

B.  Objectives 

The  objective  stated  in  the  contract  is  "to  investigate  the  prob¬ 
lems  associated  with  the  design  and  analysis  of  foundations  for  proteotive  struc¬ 
tures  which  are  subjected  to  dynanic  loads  from  nuclear  blasts".  This  general 
objective  really  indicates  little  regarding  the  technical  direction  of  the  program 
The  more  specific  goals  of  the  present  program  are  limited  to  consideration  of 
spread  footings  (this  is  contrasted  to  foundations  in  general).  At  least 
from  a  qiialitative  point'Of-view,  it  was  postulated  in  the  original  program 
that  behavior  of  all  foundations  could  be  explained  by  an  understanding  of 
the  behavior  of  spread  footings  and  pile  foundations  with  other  foundations 

^  McKee,  K.E. ,  Design  and  Analysis  of  Foundations  for  Protective  Structures 
AFSWC-TR-59->6,  15  Oct.  1959. 
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being  considered  as  some  combination  of  these  tvo.  Attention  in  this  program 
has  been  devoted  to  two  main  problem  areas  which  have  been  investigated  con¬ 
currently*  studies  of  dynamic  behavior  footings  have  been  concerned  primarily 
idth  vertical  loads  where  the  static  behavior  is  understood  relatively  well  and 
studies  of  footing  behavior  under  other  than  centrally  applied  vertical  loads 
which  to  date  have  been  carided  out  primarily  for  the  static  case. 

C,  Technical  Approach 

Ihe  approach  has  been  a  conhlnatlcm  of  analytical  and  esq^erlmental 
work.  In  a  virgin  technical  area,  it  is  desirable  that  these  two  research  tools 
proceed  together.  A  highly  developed  theoretical  approach,  regardless  of  its 
technical  sophistication  or  elegance,  must  be  able  to  withstand  the  harsh  light 
of  experimental  reality.  On  the  ot>»r  hand,  experimentation  by  itself  may 
furnish  a  mass  of  useless  data  if  theoretical  work  is  not.  available  to  provide 
a  basis  for  correlation.  A  purely  empirical  approach  is  always  expensive  and 
unless  there  Is  a  satisfactory  understanding  of  the  phenomena,  it  may  lead 
to  misinterpretations, 

D,  Diseussion  of  the  Problan 

The  behavior  of  footings  subjected  to  dynamic  loads  represents  a 
major  technical  area  in  which  little  work  has  been  done.  The  particular  objective 
of  this  project  is  somewhat  restrictive  since  consideration  is  limited  to  those 
aspects  related  to  protective  construction.  The  actual  reduction  represented  by 
this  objective  is  primarily  with  respect  to  application  rather  than  fundamental 
requirements. 

Quantitative  infornation  relating  to  the  behavior  of  snail  footings 
subjected  to  dynamic  loads  has  been  obtained  under  this  program  for  limited  soil 
types.  Additional  experimental  data  is  currently  being  obtained.  It  is  anti¬ 
cipated  that  the  analysis  of  this  data  will  provide  a  quantitative  understanding 
of  footing  behavior.  As  this  analysts  proceeds  the  direction  for  future  research 
should  become  obvious. 

Only  one  previous  experimental  study  of  the  behavior  of  footings  sub¬ 
jected  to  dynamic  loads  is  known  to  the  author.  In  19$hs  Massachusetts  Institute 
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of  Technology  conducted  limited  static  and  dynamic  footing  test8=^«  Currently 
•Wiree  other  agencies  are  testing  footings  subjected  to  dsniamlc  i'^adsi 
Naval  Civil  Engineering  Laboratoryi  Port  Huenerae,  Callfonalaj  University  of 
Illinois,  Urbana,  Illinois}  and  Waterways  Experiment  Station,  Vicksburg, 
Mississippi.  For  information  on  these  programs  the  reader  is  referred  to 
these  agencies. 

E.  Report  Organlaation 

The  !naln  body  of  this  report  is  primarily  desorlptlve  in  nature. 
Chapter  2  considers  the  theoretical  work,  while  Chapter  3  is  concejmed  with 
the  experimental  studies.  Chapter  Ut  Conoluslons,  attempts  to  generalize 
what  has  been  done  on  this  program  to  date,  and  Indicates  the  direction  the 
research  will  follow  for  the  remainder  of  the  project.  ' 

The  first  two  appendices  deal  with  plastlolty  and  limit  analyses. 
Appendix  A  presents  a  general  review  of  pertinoit  plasticity  studies  and  a 
discussion  of  their  advantages  and  llmiiiations.  Appendix  B  represents  a  first 
approach  at  applying  the  methods  of  plastlolty  to  the  dynamic  behavior  of  footings 
Appendix  C  extends  the  earlier  theoretical  atudiss  to  include  eocentrlo  loads. 
Appendix  D  considers  the  Ottawa  sand  used  in  the  experimental  program.  The 
next  two  appendices,  E  and  F,  report  on  the  two-  and  three-dlmenBlonal  footing 
tests  which  have  been  conducted  to  date*  Finally,  Appendix  0  contains  Information 
on  the  dynamic  loading  apparatus  constructed  for  this  project. 


Massachusetts  Institute  of  Technology,  The  Behavior  of  Soils  Under  Dynamic 
Loading,  3,  Pinal  Report  on  Laboratory  Studies,  AFSWP  118,  Aug  19^h.  ^ 
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Chapter  2 

TwerORKTIGAL  INVESTIGATION 


A.  General 

The  close  association  of  the  theoretical  research  with  the  laboratory^ 
Invest^ pat Inns  makes  any  clear  distinction  between  theory  and  experimentation 
nearly  Impossible  and  certainly  undes'  rable.  This  chapter  discusses  primarily 
the  theoretical  ivork  w^  t*^  reference  to  the  laboratory  studies  as  required. 

The  t'^enretlcal  work  has  proceeded  on  several  fronts.  An  attempt  will'  be  made 
to  report  propresr,  evaluate  the  suitability  of  the  approach,  and  indicate  t'^e 
direction  of  future  studies. 

3.  Plasti  city  Theory  in  Foil  Mechanics 

The  in3tabll''ty  of  the  soil  belnv;  n  footing  r.  associated  with  the 
formation  of  shear  surfaces.  Many  of  the  approaches  ur.'^cj  for  st'itii.l-' ty  'rnbl'^ms 
by  soil  mechanics  engineers  rest  on  approximate  exton:;'’.'!:).-  oj-  (rrtj'ii  c  I.  i.cai 
plasticity  solutions,  e.g,,  Prandtl's  solution  for  foof.np  failure.  The  soil 
mechanicians  from  this  basis  have  proceeded,  in  general,  with  an  engineering 
point-of-view  using  assumed  failure  surfaces  and  simplifications  to  treat  stability 
problems.  Occasional  forays  into  soil  mechanics  have  been  made  by  plastlcians 
who  have  considered  methods  of  solution  and  have  solved  specific  problems.  Neither 
the  soil  mechanics  nor  the  plasticily  theories  provide  completely  satisfactory 
solutions  neither  approach  is  able  to  accurately  predict  the  behavior  of  an 
actual  footing. 

Attention  during  the  original  program^  ras  limited  to  modifications 
and  extensions  of  more-or-less  standard  soil  mechanics  approaches.  Based  on  tbe 
original  research,  it  was  thought  that  the  more  esoteric  plasticity  approaches 
could  be  justified.  The  experimental  phases  of  the  original  program  raised  serious 
doubts  regarding  the  suitability  of  " standard' soil  mechanics  solutions  to  the 

^  McKee.  K.  E.,  OprcltV 
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the  prediction  of  the  baharior  bf  footings  subjected  to  dynamic  loads,  Ihe 
plasticity  approach,  although  complex  for  practical  applications,  gave  hope 
of  providing  a  batter  understanding  of  the  phenwnena  involved,  Appendiees  A 
and  B  are  the  result  of  these  plasticity  considerations. 


Appendix  A  contilns  a  general  discussion  of  the  plasticity  and  limit 
analysis  theory  vith  particular  attention  gl^n  to  tiie  potential  application 
to  soil  mechanics.  A  number  bf  papers  relating  the  two  subjects  are  considered 
to  demonstrate  prevdoua  vfork  in  this  area.  In  preparing  this  appendix,  an  attempt 
was  made  to  provide  an  Introduction  to  the  theory  of  plasticity  for  those  not 
well  versed  in  this  field.  Most  of  this  appendix  is  eoneemed  with  static 
problems,  with  the  concept  of  dynamic  solutions  being  introduced  only  in 
general  terms. 

Appanilx  Bis  specifically  concerned  with  an  investigation  of  a 
theoretical  approach,  based  on  the  plasticity  theory,  to  the  behavior  of  a  footing 
subjected  to  (^y^.?.m^c  loads.  The  purpose  of  the  research  was  to  consider  first, 
if  such  a  solution  were  possible,  and,  second,  to  set  up  a  method  of  approach, 

The  possibility  of  s  theoretical  approach  was  established  under  certain  restrictive 
assumptions  and  an  example  was  presented.  Greater  generalization  to  extend  the 
solution  to  other  problems  seems  to  be  possible  although  the  actual  labor  of 
obtaining  such  solution:  is  expected  to  Increase  substantially,  IVie  oveir-all 
results  of  this  appendix  are  most  Interesting  since  the  possibility  of  t)i4s  type 
cf  theoretical  approach  has  been  established.  From  the  point-of-vlew  of  long- 
range  research  into  the  behavior  of  footings  under  dynamic  loads,  this  typs  of 
approach  may  provide  the  final  answer  and  additional  studies  in  this  area  should 
be  considered. 


C,  Dynamic  Response 


The  behavior  of  footings  subjected  to  dynamic . loads  was  analyzed 
in  the  original  report  by  what  m.ight  be  termed  an  "engineering  approach", 

Th'' ?  aprroach  barpd  pn  an  extension  to  tine--dependpnt  loads  of  Andersen's^' 


Andersen, 
New  Tcr'Xo 


P, :  Ptubstructure  Analysis  and  Design,  The  Donald  Press,  New  YorV, 

prT'7. 


ARMOUR  RESEARCH  FOUNDATION  OF  LiNOIS  INSTITUTE  OF  TECHNOtOGT 


ARF  Project  No.  8193-7 
Interim  Report 


work  relating  to  one-sided  footing  failure.  The  major  assumptions  introduced 
for  this  approach  arei 


1,  The  failure  surface  under  (^mamlc  loads  will  be  ■Wie  surface 
determined  by  application  of  the  initial  value  of  the  ovei> 
pressure  as  static  pressure, 

2,  The  resl stance  offered  to  movement  ia  rigid  plastic  in  form, 
i,e,,  settlement  and  soil  conpresslbility  are  not  considered, 

3,  The  maximum  plastic  resistance  equeils  the  static  resistance 
determined  analytically,  assuming  tiie  failure  surfaces  considered 
are  the  a4»me  as  in  item  1  above, 

li.  The  behavior  of  the  soil  Is  governed  by  the  parameters  7,0, 

and  c,  where  0  and  c  may  themselves  be  functions  of  many  para¬ 
meters  relating  to  the  soil  and  the  condition  of  loading. 


By  considering  tho  notion  of  the  failure  mass  of  soil,  a  differential  equation 
was  established  for  the  dyToamle  behavior.  Solutions  of  this  equation  allow 
prediction  of  the  dlsplacerients,  investigation  of  inertia  effects,  etc. 

In  the  original  program,  the  one-sided  fallurs  pattern  CAndersen’s) 
was  specially  selected  because  it  was  compatible  with  the  dynamio  analysis.,^ 

This  canpares  with  other  fallurs  assumptions,  such  as  Terzaghi’s  modification 
of  Prandtl's  solution,  for  which  the  failure  mode  is  oompatible  only  for 
incipient  motion.  Two-sided  failure  patterns  generally  are  not  suitable  for 
gross  motion.  Since  two- sided  failures  do  occur,  an  attempt  has  bean  made 
to  consider  possible  two-sided  failure  patterns,  The  approaches  considered 
are  discussed  in  the  following  paragraphs, 

(1)  The  first  attempt  to  arrive  at  a  suitable  approach  was  based  on  the 
following  assumptions  (refer  to  Fig,  1  for  nomenclature)! 

1,  A  wedge  of  soil  under  the  footing  Is  assumed  to  remain  rigid 
and  to  move  with  the  footing  (defined  by  sides  making  an  angle 

of  0  with  the  horizontal  surface  for  footingB  with  a  rough  base), 

2,  The  failure  surface  is  a  segment  of  a  circle  from  the  apex  of 
the  wedge  to  the  surface. 

3,  The  center  of  the  failure  circle  is  located  arbitrarily  at 

a  distance  y  above  the  surface  and  at  a  horizontal  distance 
X  from  the  centerline  of  the  footing, 

h.  The  failure  circle  associated  with  the  lowest  failure  load, 

P  ,  is  assumed  to  be  the  most  nearly  correct  circle. 

5,  The  failure  Is  symmetrical  about  the  centerline  of  the  footing, 

6.  The  full  cohesive  stress  is  developed  in  the  soil  along  the 
failure  surface. 
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7.  The  frictional  fore©  along  the  failure  arc  is  governed  by  the  cohesion 
and  the  velght  above  each  point. 

These  assumptions  can  be  Justified  for  static  conditions  in  at  least  ttro  vayai 
first,  the  resulting  failure  surfaces  are  similar  to  those  postulated  by 
Terzaghi  for  a  rough  based  footing^,  and,  second,  the  failure  surfaces  are 
similar  to  those  obtained  from  footings  on  Ottawa  sand^,  Ihese  Justifications 
are  by  no  means  introduced  to  prove  that  the  assumed  failure  pattern  is  correct. 

It  is, however,  felt  that  the  aeaumed  failure  pattern  will  be  a  reasonable 
representation  for  the  observed  pattern^.  Also,  since  no  other  assumption, 
including  Teraaghl'  si/,  is  oonpletely  In  agreawnt  with  the  observed  patterns, 
one  might  expect  the  asssumptlons  made  here  to  provide  essentially  the  same 
bearing  capacity  as  Terzaghi 's  solution  for  static  loads, 

(2)  The  second  approach  represents  a  modification  of  Hill's  plasticity 
solution  (see  Appendix  A).  For  this  approach,  no  wsdge  is  assumed  under  the 
footing,  and  the  failure  pattern  is  as  indicated  on  Fig,  2»  If  one  considers 
this  failure  mode  as  two  mirror-image  footings  exhibiting  one-sided  failure, 

the  approach  used  In  the  original  program  for  one-sided  failure  (based  on 
Andersen's^  fsilure  surface)  is  directly  applicable,  If  the  footing  Is  oon- 
sidered  as  half  the  width,  the  one-sided  failure  load  can  be  doubled  to  ob¬ 
tain  the  bearing , capacity  for  the  original,  footing,  The  nature  of  this  approach 
is  so  similar  to  that  used  earlier  for  one-sided  fad. lure  that  detailed  con¬ 
sideration  is  not  required.  Plots  of  P  versus  B  or  r  versus  B  could  be 
easily  prepared  for  various  soil  properti.es  as  was  done  earlier.'^  In  considering 
this  approach  it  should  be  kept  in  mind  that  in  the  experimental  phases  of  this 
program  wedges  are  observed  under  the  footing. 

(3)  The  plasticity  solvUtion  for  the  punch  on  the  half  space  (Prandtl's 
solution)  can  also  be  used  directly.  Figure  B,3  (Appendix  B)  shows  the  velocity 
field  associated  with  the  Prandtl  solution.  Such  a  failure  mode  could  be  used 
directly  as  a  basis  for  predicting  symmetrical  failure, 

Terzaghi,  K,  and  Peck,  R.  B.,  Soil  Mechanics  In  Fhgineering  Practice, 

John  Wiley  and  Sons,  Inc.,  New  York,  19 hd,  p.  1^(J, 

McKee,  K,  S.,  Op.  Cit,,  Appendix  C 

1/ 


V 


Andersen,  P.,  Op,  Pit.,  p,  81 


McKee,  K,  E,,  Op,  Pit.,  Appendix  S 
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These  three  approiehes  should  be  considered  only  as  examples  of 
the  many  failUTe  .j'^’tterns  which  might  be  postsilated.  AllAough  soro  con¬ 
sideration  has  been  given  to  each  of  the  three  approaches,  they  are  limited 
by  their  common  aastmptlons.  Detailed  research  in  this  area  will  be  oorrfucted 
when,  and  if,  it  is  justified  by  the  eaqperlmental  results.  For  the  present, 
the  aim  was  to  consider  some  possible  method  of  treating  two-sided  failure  with 
the  sane  general  type  of  engineering  approach  used  earlier  for  ono-slded  failure. 

It  may  be  well  to  consider  some  of  the  limitations  implied  in  the 
above  approaches,  with  emphasis  on  the  observed  experimental  results.  The 
failure  mode  under  dynamic  loads  may  be  the  same  as  for  etatio  loads,  i.e., 
a  clearly  defined  failure  surface,  but  the  experimental  results  for  Ottawa  sand 
indicate  that  such  a  failure  mode  may  form  for  dynamic  loads  only  when  dis- 
placaments  are  much  greater  than  would  be  associated  with  their  fonnatlon  under 
static  loads.  This  indicates  that  the  failure  mode  may  differ  from  the  Idealized 
assumptions  discussed  above  through  at  least  part  of  the  displacement.'  The  actual 
resistance-displacement  is  not  rigid-plastic  for  static  loads,  and  tJMire  is  no 
reason  to  assume  that  It  would  be  under  dynamic  conditions.  The  other  elements 
of  the  assumptions  similarly  must  be  subjected  to  critical  evaluation  before 
they  can  be  accepted.  Suffice  it  to  say,  that  there  are  many  questions  which 
only  the  evaluation  of  the  experimental  results  can  answer, 

D.  Non-Vertical  Loads 

In  the  original  program,  attention  was  devoted  primarily  to  footings 
subjected  to  vertical  forces.  This  is  an  extremely  common  type  of  loading, 
and  for  this  reason  it  can  be  expected  to  be  important,  Thera  are,  however, 
other  types  of  loads  which  do  occur  alona  or  in  combination  with  the  vertical 
loadsi  (1)  overtum:^ng  moments,  (2)  inclined  loads,  and  (3)  torsional  moments. 

The  last  of  these  is  of  little  practical  importance  and  will  not  be  considered 
in  this  report. 

Overturning  moments  can  be  treated  by  considering  eccentric  loads, 
i.e.,  combinations  of  direct  vertical  load  and  overt'oTning  moments.  3y  intro¬ 
ducing  the  same  assumptions  used  in  the  original  analysis  for  one-sided  failure, 
this  problem  can  be  analyzed.  ■  The  analysis  for  this  class  of  loading  is  given 
in  detail  in  Appendix  C,  Since  the  treatment  there  is  complete,  suffice  it 
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to  say  -ttiat  tiie  graphical  results  presented  in  Appendix  E  of  the  final 
reporw  can  be  used  with  only  superficial  modifioationa.  It  should  be 
noted  that  for  this  ourpose  the  one-sided  failure  mo*  represents  the 
actual  behavior. 

The  influence  of  inclined  loads  has  been  studied  during  the 
present  program,  but  to  date  the  results  are  not  satisfactory.  It  is  anti¬ 
cipated  that  further  consideration  will  result  in  more  suitable  analytie 
aoproaehes  and  furthermore  should  allow  for  a  suitable  form  of  presentation. 
At  the  time  of  this  mrlting,  the  need  for  a  method  of  treating  footings  sub¬ 
jected  to  inclined  loads  is  pointed  out,  and  the  fact  that  work  has  been  and 
will  be  done  in  this  respect  ia  reported. 


IT - 

McKee,  K.  b., 
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Chapter  3 

SCPEHIMENTAL  INVESTIGATION 


A.  General 

The  scope  of  this  project  Includes  latoratory  e^erlaentation 
to  verify  and  support  the  theoretical  research.  This  chapter  discusses 
the  planned  experimental  studies  and  reports  on  the  results  obtained  thus  far. 
.*?lnca  the  laboratory  work  was  initiated  during  June,  the  ej^eriinsntal  results 
are,  of  necessity,  limited,  and  In  some  cases  only  preliminary  results  are 
available.  Tn  the  following  sections,  the  experimental  studiea  will  be  con¬ 
sidered  in  three  phaaoai  material  properties,  twoHilmenslonal  experiments 
In  the  glass  box,  and  three-dimensional  experiments  In  the  sand  box.  The 
planning  for  each  of  these  phases  will  be  diBCuaaad  with  detailed  infonnatlon 
reported  in  the  appendices, 

S,  Material  Properties 

Material  property  studies  fall  into  two  categorleat  standard 
classification  tests  required  as  part  of  the  over-all  experimental  pxxjgram 
and  special  research  intended  to  increase  the  available  kncwledge  regarding 
the  properties  of  dense  sand.  The  first  Is  an  adjunct  to  the  experimental 
resetu'ch  and  requires  no  special  justification.  Suffice  It  to  say,  tiiat 
standard  soil  classification  tests  will  be  conducted,  supplemented  by  addi¬ 
tional  special  tests  as  may  prove  desirable.  The  second  category  consists  of 
special  studies  ijitended  to  provide  detailed  soil  properties.  Appendix  D  reports 
on  laboratory  work  dealing  with  Ottawa  sand.  The  angle  of  internal  friction, 
found  in  the  triaxlal  tests  of  the  original  progaram  was  questioned  since  it 
differed  significantly  from  that  computed  from  bearing  formulas  based  on  the 
footing  tests  conducted  in  the  laboratory.  Improved  experimental  procedures 
have  provided  data  regarding  the  angle  of  internal  friction  —  data  which  appears 
to  correlate  with  the  results  of  the  footing  tests.  Also  included  in  Append-' x  E 
are  the  results  of  a  series  of  experiments  dealing  wj  th  restraints  at  the  inter¬ 
face  between  sand  and  glass  or  sand  and  steel,  These  studies  are  intended  to 
pro-vlde  relative  and  absolute  information  for  use  in  the  glass  box  (two- 
dimensional  experiments). 
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For  the  initial  laboratory  studies  involving  cohesive  soils,  it  was 
decided  to  use  one  of  the  manufactured  soils  that  have  been  investigated  under 
ARF  sponsorship^^.  By  combining  natural  soil  constituents  under  controlled 
conditions  it  has  been  found  possible  to  produce  soils  having  consistent 
prcperties  suitable  for  experimental  studies, 

C.  Two-»Dimens1  onal  Studies 

A  glass-sided  container  has  been  used  to  slnulata  two-dlnansional 
conditions.  Static  and  dynamic  tests  have  bean  conducted  on  loose  and  dense 
Ottawa  sand,  and  the  manufactured  cohesive  material.  Static  loads  were 
applied  through  a  hydraulic  jack  or  by  means  of  a  controlled  strain  rate 
apparatus.  Dynamic  loads  vmr©  applied  using  the  apparatus  developed  on  this 
program  (see  Appendix  0).  Complete  records  were  made,  and  rart-’x  pV'Otcgranhc 
ware  taken  of  t''?©  soll-foof’ np  interaction  under  dynamic  loads.  Detg^l8  of 
these  experiments  are  contained  in  Appendix  B, 

D ,  Three-Dimensional  Studies 

A  box  having  plan  dimensions  of  )i  ft  by  la  ft  and  a  depth  of 
2  ft  was  used  for  three-dimensional  footing  tests.  The  only  material  utilized 
to  date  has  been  Ottawa  sand.  This  same  material  war  used  in  the  original  program, 
thus,  the  results  of  the  r*atic  experiments  conducted  in  the  original  program  are 
applicable.  Details  of  the  three-dimensional  dynamic,  experiments  are  re¬ 
ported  in  Appendix  FI 


^  Sellg,  E,  T.,  and  Rowe,  R.  D.,  Artificial  Foils,  Armour  Research  Foundation, 
K920,  September  1960, 
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Chapter  li 
COITOLIBTOIB 


This  Interim  report  has  attempted  to  cover  the  research  conducted 

during  t>’R  first  half  of  the  contract  period.  In  peneral,  the  reported  work 
is  in  profrress,  and  as  a  result  complete  data  and  generallaations  basad  on 
the  research  are  lacking.  Thus,  the  conclusions  and  recommendations  expressed 
below  are  tentative  and  sub,1oct  to  later  revision. 

The  partnership  of  theoretical  and  experimental  researoh  applied 
-  on  this  project  has  proved  most  effective.  Parallel  devolopnent  in  ttiese  two 
directions  has  maximized  the  information  which  has  been  obtained  within  -010 
scope  of  the  research.  The  availability  of  experimental  results  has  lead  to 
the  elimination  of  several  'premising'  theoretical  approaches,  whils  at  the 
same  time  observationa  have  provided  the  concepts  for  other  analytical  techniques. 

The  studies  to  date  have  provided  much  information  regarding  footing 
behavior.  The  studies  relating  to  behavior  under  static  loads  verify,  in 
general,  the  classical  shear  surface  approaches.  One-sided  and  eymmetrlcal 
failure  surfaces  have  been  observed.  The  bearing  capacities  have  been  correlated 
vrtth  those  determined  from  tho  applicablo  theories.  Earlier  incongniitieB, 
between  the  soil  properties  determined  from  triaxlal  testa  and  those  required 
to  correlate  the  theories  and  experiments,  have  been  resolved.  All  Indications 
are  that  txiaxial  or  direct  shear  tests  conducted  on  the  same  soil  and  at  the 
same  density  vjould  result  in  soil  properties  satisfying  the  theoretical  requirements 

The  original  project^  considered  only  one  type  of  dynamic  load  — 
the  dropped  weight.  These  experiments  were  carried  out  to  establish  the 
potential  of  conducting  dynamic  experiments  on  small  footings,  but  due  to 
project  limitations,  recording . was  limited  to  visual  and  photographic  observa¬ 
tions  before  sind  after  the  load  application.  From  these  observations  of 
for!tinf:s  in  both  the  sand  box  and  the  glass  box,  there  was  no  indication 
of  the  formation  of  a  shear  surface  under  a  single  load  application.  (The 
second  and  third  drops  on  the  same  footing  were  sufficient  to  cause  shear 
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failure.)  This  failure  mode,  i.e.,  settlement  without  observable  shear  surfaces 
under  dynamic  loading  caused  by  a  single  dropping  of  a  weight,  would  bo  quite 
diffewmt  from  that  observed  at  sane  defleotlon  for  a  static  load.  This 
'c!if;p'srent'  behavior  coild  be  at  least  intuitively  explained,  since  the  dropped 
weight  represents  a  relatively  extreme  dynamic  load  —  es^ntlally  an  initial 
impulse.  On  the  other  hand,  it  could  be  argued  that  a  drop  ftcsn  a  sufficient 
height  would  cause  shear  surfaces  to  fern  and  hence  that  the  behavior  would  bo 
similar  to  that  for  static  loads. 

The  development  of  the  dynani c  loading  apparatus  on  the  current 
project  has  allowed  the  application  of  less  extreme  ^pe  of  dynamic  loads.  The 
-results  of  such  testa  Indleato  that  shear  surfaces  would  form  although  the 
associated  displacements  exceed  those  for  static  load.  This  Implies  that  the 
more  extreme  the  type  of  dynamic  load,  the  greater  the  difference  in  behavior, 
(in  the  case  of  a  linearly  decreasing  load-time  history,  extremity  would  bo 
associated  \rtth  high  peak  loads  and  s‘-ort  durations).  It  might  further  be 
speculated  that  the  difference  In  behavior  is  associated  with  the  displacement 
at  wh<  eh  failure  e^lr faces  fbrm. 
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Appendix  A 

APPLICATION  OF  PIASTICITT  AND  LIMIT  ANALYSIS 
TO  FOOTINO  TE3ION 

by  R.  L.  Chiapatta 
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Aooonriix  A 


A^PI.Tr.ATTQ?!  OF  T^LASnCTTy  AND  LT^rrT  MlAUSir 

TO  r'oonx  D’^sicr 


by  R.  L.  Chiapetta 


A . 1  introduction 

Thr  objoctivo  of  t^^ia  appondix  is  to  consider  the  rolo  of  the 
r'-oory  of  olastic^  ty  and  lin^t  analysis  as  tools  for  invnsf' rratl  np;  the 
bohavior  of  footlnpe.  The  behavior  of  footinpe  subjected  to  static  loads 
v;ill  bo  ri-trcupsod  :)nrticular  attention  to  the  dlfficultios  involved 

in  obtaininr  a  solution  and  aomo  of  the  limitations  of  tho  methods.  This 
ulll  be  follov;Qd  by  a  proFontation  of  results  of  past  studies  of  tho  static 
oroblom.  Consideration  of  the  dynamic  problem  la  contained  in  Appendix  B, 
^plld  tly,  the  discussion  will  be  mainly  restricted  to  consideration  of  con¬ 
tinuous  footlnpe  subject  to  vertical  loads.  It  is  assumed  that  the 
footing  la  rifH  d  relative  to  t^o  soil,  so  that  the  footinp  may  be  considered 
afj  a  ripld  punch. 

The  difficulty  in  determining  the  behavior  of  a  footing  subjected 
to  either  static  or  dynamic  loads  stems,  in  part,  from  the  complex  mechanical 
properties  of  real  soils.  For  simplicity,  the  soil  properties  are  idealized  to 
perm-*  t  use  of  mathematical  theories,  such  as  the  theories  of  elasticity  and 
plasticity.  Both  of  these  theories  may  be  applied  to  the  footing  response 
problem,  ^Infortunately,  the  t^’eor:/’  of  elnstl  city,  at  times,  has  been 
employed  where  its  validity  is  questionable.  This  is  due  in  part  to  the 
absence  of  an  adequate  formulation  of  a  rule  for  determining  the  validity 
of  t^'e  theory,  Th^  follov;ing  "rule  of  thumb"  for  the  applicability  of 
elastic  theory  in  soil  mechanics  problems  has  been  suggested  (l)i/.  If  the 
factor  of  safety  of  a  mass  of  soil  with  respect  to  failure  by  plastic  flow 
(continents  Hofcrr.ation  at  a  constant  st-cto  of  stress)  exceeds  a  value 
of  about  3,  tbe  vertical  stresses  in  the  soil  can  be  estimated  by  the  theory 
of  einrticHy. 

Numbers  in  parenthesis  refer  to  references  listed  in  Section  A, 7, 
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A, 2  Plasticity  Theory 

In  this  appendlXi  failiare  ia  defined  such  that  In  some  portion  of 
the  soil  the  stress  level  is  hlRh  enough  to  extend  into  the  "plastle"  range# 

It  is  further  assumed  that  the  stresses  at  a  point  are  related  according  to  an 
accepted  failure  theory.  For  failure,  it  is  expected  that,  in  general,  three 
types  of  zones  in  the  soil  will  have  to  be  considered:  elastio  zones  where 
the  stresses  can  be  conputad  by  elaotie  theory,  plastic  zones  whore  the  state 
of  stress  may  be  determined  by  using  equations  of  plasticity,  and  transition 
zones  of  Intermediate  states  of  stresa.  The  existence  of  transition  zones 
makes  the  problem  of  computing  the  stresses  eflctranely  complicated.  In  order 
to  simplify  the  analysis,  the  existence  of  transition  zones  la  sometimes 
disregarded,  and  the  soil  is  referred  to  as  an  alastio-plastla  material. 
Although  swie  work  has  been  conducted  on  work-hardening^  theories  for 
soil  mechanics  (2),  the  theories  have  not  been  developed  to  the  point  of 
general  acceptability.  For  this  reason  it  is  usually  assumed  that  the 
material  is  perfectly  plastic.^ 

The  field  oquationa  of  plasticity  consist  of  equilibrium  equations, 
atraln-dieplacement  aquations,  stress-strain  relations,  and  the  yield  funotlon. 
An  exact  solution  to  a  plasticity  problan  consists  of  a  stress  field  and  a 
corresponding  strain  or  displacement  ^eld  related  by  the  stress-strain 
relations.  The  stresses  must  satisfy  the  equilibrium  equations,  the  yield 
function,  and  boundary  conditions  on  the  stresses.  The  compatible  dis¬ 
placements  must  satisfy  specified  displacement  boundary  conditions. 

The  equilibrium  and  strain-displacement  equations  are  of  the  same 
form  as  in  the  theory  of  elasticity.  However,  in  plasticity  the  total  strains 
appearing  in  the  strain-displacement  relations  are  composed  of.  the  sum  of 
the  elastio  and  plastic  strains  (3).  The  elastio  strains  are  related  to  the 
stresses  by  the  stress-strain  relations  of  linear  elasticliy.  The  usual 
assumptions  of  perfect  plasticity  require  ihe  stress-plastic  strain  rate 


A  work-hardening  material  is  one  for  which  increased  stress  accompanies 
Increased  strain  in  the  plastic  range. 


2/ 


A  perfectly  plastic  material  ie  one  for  which  the  strain  will  Increase 
under  constant  stress. 
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relation  to  be  of  the  form 


1,1 


where  the 


^  3<rij  ■ 

are  the  comoononts  of  the  nlastic  strain  rates, 


(A.l) 


are  tne  comoononts  oi  tne  niastic  strain  rates,  arc 

the  correspondinE  comnoncr.ts  of  stress,  F  is  the  yield  function,  and  ^  i: 
a  positive  factor  of  oronortionality  which  my  be  a  function  of  the  space 
coordinates  (U). 


In  tv/o- dimensional  soil  mechanics.  Coulomb's  formula  (1)  for  the 
yield  function  is  usually  assumed.  This  function  reduces  to  the  Tresca  or 
von  tiises  function  v/hen  the  an/tle  of  internal  friction  is  equal  to  aero. 

’rfith  reference  to  three  dimensions,  the  Coulomb  failure  law  for  an  ideal  „  - 

cohesive  soil  has  been  used  to  obtain  the  yield  surface  reoresenting  the 
yield  function  (5).  Modified  Tresca  and  von  Jiises  yield  functions,  which  are 
denendent  on  mean  normal  stress,  have  been  considered  as  proner  generalizations 
of  the  Coulomb  rule  to  three  dimensions  which  reduce  to  the  Coulomb  law  for 
two  dimensions  (6),  An  imortant  feature  of  a  nerfectly  elastic  material, 
that  obeys  Coulomb's  failure  theory,  is  th.at  it  exhibits  dllatancy,  that  is, 
clastic  deformation  must  be  accomnanied  by  an  increase  in  volume  (fl).  It  should 
be  noted  that,  in  fact,  dilatancy  is  observed  in  many  soils,  such  as  dense 
sand  or  stiff  claor. 


For  the  problem  cf  a  rigid  nunch  on  a  semi-infinite  mass  of  soil, 
only  the  incioient  clastic  flow  oroblem  is  usually  considered  in  the  literature, 
so  that  the  boundary  oonJitions  at  the  undeformed  surface  are  satisfied  (7). 
Determining  the  stresses  and  velocities,  after  the  ounch  has  penetrated  a 
finite  distance,  is  more  difficult  and  would  require  a  study  of  the  successive 
phases  of  the  plastic  flow.  A  consideration  of  the  static  problem  of  footings 
Is  helpful  in  gaining  insight  into  and  is  a  necessary  prelude  to  the  stuc^  of 
the  corresponding  dynamic  problem. 

The  case  of  a  continuous  footing  reduces  to  a  plane  strain  problem. 

Ln  general,  for  an  elastic-perfectly  plastic  material  three  ranges  of  mechanical 
behavior  are  present  in  a  olane  strain  oroblemi  the  elastic  range  (below  the 
elastic  limit),  the  range  of  contained  clastic  deformation  (between  the  elastic 
limit  and  the  flow  limit),  and  the  range  of  unrestricted  clastic  flow  (beyond 
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the  flow  liraLt),  t-Iany  plsuie  strain  problams  have  been  solved  by  'Uie  theory 
of  elastleltyi  but  analytical  solution  for  the  entire  range  of  contained 
olastlo  deform tliwi  has  generally  been  possible  only  when  the  slvipe  of  the 
elastlc'plastlc  boundary  was  known  a  priori  from  symuBtry  oondltlonst  As  a 
rule,  the  shape  of  the  alastio-plastlo  boundary  is  not  known  beforehand  and 
general  methods  of  8olTd.ng  such  problems  have  not  been  dovslopad  (3)i  An 
additional  complication  is  that  the  stress  boundary  conditions  are  not 
generally  sufficient  to  make  the  problem  statically  determinate. 

Because  of  these  difficulties,  contained  plastic  deformations  can 
seldom  be  treated  by  exact  analy^s,  and  it  is  necessary  to  consider  numerical 
methods  such  as  have  been  successfully  used  in  the  case  of  elastlo-plastio 
torsion.  However,  locating  the  elastic -plastic  boundary  by  numerical 
techniques  is  much  more  difficult  in  plane  strain  than  in  torsion  (3). 

Although  a  stress  function  describes  the  stress  distribution  in  both  cases, 
the  stress  components  are  given  by  the  first  derivatives  in  the  case 
of  torsion,  but  by  the  second  derivatives  in  the  case  of  plane  strain. 

In  both  cases,  the  elastic-plastic  boundary  is  located  by  the  condition 
that  dB.acontinultles  of  stress  cannot  exist  across  this  boundary.  In  the 
case  of  torsion  this  means  that  the  stress  function  and  its  first  derivative 
should  be  continuous  across  the  olastic-plastlc  boundary.  However,  for 
plane  strain,  the  stress  function  and  its  first  and  second  derivatives 
must  be  continuous.  To  fulfill  the  condition  of  continuity  at  the  elastic- 
plastic  boundary,  a  much  finer  mesh  must  be  used  in  problems  of  plane 
strain  than  is  necessary  in  torsion  oroblems.  This  can  be  expected  to  oom- 
ollcata  the  numerical  work  considerably. 

The  final  state  of  unrestricted  plastic  flow  usually  cannot  be 
treated  in  plane  strain  problems  without  a  full  analysis  of  the  precetS-ng 
states  of  ccautained  plastic  deformation.  The  analysis  of  the  contained 
plastic  deformation  would  be  necessary  to  locate  the  boundary  between  the 
elastic  and  plastic  regions.  In  the  problem  of  torsion  of  bars,  it  is 
possible  to  treat  the  final  state  of  unrestricted  plastic  flow  without 
reference  to  the  preceding  elastic-plastic  states  since  such  flow  can  occur 
only  when  the  entire  bar  has  become  plastic  (3). 
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A. 3  limit  Theories 


If  one  is  interested  in  the  loads  necessary  to  cause  unrestricted 
olastic  flow  and  not  particularly  interested  in  the  o:act  stress  distribution 
at  failure,  the  thaorenia  of  limit  desipn  (9)  can  bo  a  prsat  aid.  It  should 
be  noted  that  t^o  theorems  are  valid  for  three- dimonsional  problems  as  well 
as  for  those  of  olane  strain.  The  theorems  of  limit  design  are  concerned 
with  bodies  of  perfectly  plastic  material  subject  to  arbitrary  histories 
of  loading  that  are  completely  specified.  It  is  assumed  that  the  boundarj’ 
conditions  are  of  the  stress  type,  that  is,  at  every  point  of  the  body 
each  component  of  the  surface  traction  is  specified  except  where  the 
corresponding  component  of  dlsplaceiriant  is  prearlbed  to  be  zero. 

Before  quoting  the  theorems,  it  may  be  wall  to  define  some  concepts 
and  tarns  involved  in  tho  statement  of  the  thoorernsJ 

A  perfectly  oLaatic  material  is  comelotoly  characterized  by  its 
yield  function  F  .which,  for  a  homogeneous  material,  is  only  explicitly 
denendent  on  the  nine  stress  components.  States  of  stress  within  tho 
elastic  range,  0<F<1  ,  are  called  safe.  Plastic  flow  can  occur  only 
wider  states  of  stress  for  which  F  •  1  .  States  of  stresses  for  which 
F  >1  are  not  possible  in  a  perfectly  olastic  naterial, 

The  term  collapse  in  the  theorems  refers  to  conditions  for  which 
olastic  flow  would  occur  under  constant  loads  if  the  accompanying  change 
in  the  geometry  of  the  structure  or  body  were  disregarded.  For  this 
type  of  collapse,  the  equilibrium  conditions  can  be  set  up  for  the  un- 
defomed  body. 

Statically  admissible  stress  fields  are  defined  as  states  of  stre® 
for  which  the  components  of  the  stress  tensor  are  continuous  functions 
satisfying  the  conditions  of  static  equilibrium  throughout  the  botfy  and 
on  those  portions  of  the  boundary  where  the  components  of  the  surface 
tractions  are  given.  The  preceding  definition  may  be  generalized  to 
include  stress  fields  with  a  finite  nunter  of  surfaces  of  discontinuity. 

In  these  cases,  the  stresses  must  satisfy  the  conditions  of  equilibrium 
on  either  side  of  such  a  surface  and  the  components  of  the  surface  tractions 
must  be  continuous  across  the  discontinuity  surface. 

A  kinematically  admissible  velocity  field  must  satisfy  two  basic 
requirements.  'First,  the  velocity  component  must  vanish  on  those  portions 
of  the  surface  where  the  corresponding  coiroonents  of  the  surface  tractions 
are  not  orescribed,  and  secondly,  the  rate  at  which  the  actual  surface 
tractions  and  body  forces  do  work,  based  on  this  velocity  field,  must 
equal  or  exceed  the  rate  of  Internal  dissipation  of  energy  computed  from 
the  strain  rates  treated  as  ourely  plastic  strain  rates. 
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Subject  to  the  above  asaunrotions  and  definitlona  the  followinp 
two  theorems  are  «lj,di  I.  If  a  safe  statically  admissible  state  of  stress 
can  be  fo\»id  at  each  stage  of  loading,  collapse  will  not  occur  under  the  piven 
loading  history.  II.  If  a  kinematically  admissible  collapse  state  can  be 
found  at  any  stage  of  loading,  collaose  must  be  impending  or  must  have  taken 
place  previously.  These  theorems  can  be  used  to  provide  an  upper  and  lower 
bound  on  conditions  at  collapse.  The  bounds  so  establishod  may  bo  suffi¬ 
ciently  close  for  engineering  usage.  In  sone  cases,  as  will  be  mentioned 
later,  the  uoper  and  lower  bounds  can  be  made  tc  cdlneide,  thus  the  exact 
collapse  condition  is  determined. 

In  the  above  definition  of  a  klnenatlcally  admissible  velocity 
field,  no  mention  Is  made  of  continuity.  The  concept  of  a  disoontinuity 
in  the  velocity  field  often  is  very  useful  in  applying  the  limit  thaorems, 

A  velocity  disoontinuity  is  slranly  an  idealization  of  a  continuous  distri¬ 
bution  in  which  the  velocity  changes  very  rapidly  across  a  thin  transi'U.on 
layer.  If  the  yield  function  decenda  on  the  mean  normal  stress  (as  it  does  in 
Coulomb's  failure  theory),  a  discontinuity  in  tangential  velocity  must  be 
accomnanied  by  a  senaration  or  discontinuity  in  normal  velocity.  In  such  a 
case,  the  actual  transition  layer  must  have  aporeciable  thickness,  bat  the 
conceot  of  a  discontinuity  surface  still  may  be  useful  for  ourposee  of. 
calculation. 

The  above  limit  theorems  are  valid  in  the  oresence  of  a  transition 
layer  and  are  therefore  valid  in  the  limit  as  the  thickness  of  the  transition 
layer  aporoaches  zero,  provided  that  the  rate  of  dissipation  of  energy  in 
the  transition  layer  aoproaches  a  finite  limit  (9). 

A  very  simple  tyoe  of  velocity  field  that  is  often  used  for 
establishing  a  kinematically  admissible  field  is  a  rigid  body  sliding  between 
two  portions  of  the  soil  mass.  This  of  course  implies  a  disoontinuity  in 
velocity  at  the  interface  of  relative  motion.  This  type  of  disoontinuity 
will  be  called  a  "slide"  discontinuity.  As  was  pre-ziously  mentioned,  for 
a  material  whose  yield  function  depends  on  the  mean  normal  stress,  a  dis¬ 
continuity  in  tangential  velocity  requires  an  accompanying  normal  or  separation 
velocity.  For  a  slide  discontinuity  in  a  Coalonto  soil,  the  angle  between  the 
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resultant  velocity  vector  and  the  tangent  to  disoontlnulty  surface  is  equal 
•to  the  angle  of  Internal  friction.  This  condition  restricts  the  types  of 
admissible  ''slide"  discontinuity  surfaces  to  be  either  a  plane  (translation) 
or  logarithmic  spiral  (rotation).  For  a  material  whose  yield  function  is 
Independent  of  the  mean  normal  stress,  such  as  the  von  Mises  yield  fuiwtion, 
the  permissible  slide  discontinuity  surfaces  are  nlanes  and  sepnenta  of  ciroles. 
An  example  of  the  use  of  a  plane  velocity  dlsa)ntlnulty  surface  will  be  given 
Inter,  in  tHo  problem  of  determining  an  upper  bound  on  the  critical  height 
of  an  unsupported  vertical  bank  of  soil. 

Discontinuous  states  of  stress  as  well  as  discontinuous  velocity 
fields  are  permissible  in  the  application  of  the  limit  theorems.  Discontinuous 
stresa  fields  are  of  special  value  when  used  to  obtain  lower  bourrfs  on  the 
quantities  at  collapse,  even  though  the  stress  fields  may  be  without  physical 
significance  (lO).  Equilibrium  considerations  require  that  normal  and  shear 
stresa  be  continuous  across  a  line  of  stress  discontinuity,  but  the  tangential 
stress  may  be  discontinuous.  It  can  be  shown  that  a  failure  line  cannot  be 
a  line  of  discontinuity  in  the  stresses  and  since  a  discontinuity  in  the 
velocity  can  only  occur  across  a  failure  line  (10),  It  follows  that  the 
velocity  field  must  be  continuous  across  a  line  of  stress  discontinuity. 

It  has  been  sho'wn  (11)  that  the  limit  theorems  previously  stated 
for  assemblagss  of  perfectly  plastic  bodies  do  not  always  apply  when  tiiere 
is  finite  sliding  friction  at  the  common  interfaces,  ^uch  as  might  be  present 
at  the  contact  area  of  the  rigid  footing  and  the  soil  mass  In  the  classical 
two-dimensional  footing  problem.  Since  the  displacement  or  velocity  vector 
makes  an  angle  equal  to  the  angle  of  internal  friction  with  the  sliding 
surface,  there  is  a  force  acting  norral  to  the  surface,  and  work  is  dona, 
therefore,  against  the  normal  force.  It  is  this  negative  work  which  gl'ves 
rise  to  trouble  in  the  upper  bound  theorem.  The  rate  of  internal  dissipation 
cannot  be  calculated  in  all  cases  because  frictional  dissipation  is  not  uniquely 
determined  by  *:hR  flow  pattern.  It  depends  not  only  cn  relative  displacement 
rate  but  also  on  the  formal  pressure  on. the  frictional  interface,  a  quantity 
which  will  often  not  be  known.  This  type  of  difficulty  does  not  a.ffect  the 
lower  bound  theorem. 
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Theorems  have  been  developed  (11)  which  relate  the  limit  loads 
with  finite  Coulomb  friction  to  the  extreme  cases  of  zero  friction  or  of 
comolete  attachment.  The  theorems  are  stated  as  follows:  A,  Any  set  of 
loads  that  oroduco  collaose  for  the  condition  of  no  relative  motion  at  the 
Interfaces  will  oroduoe  collapse  for  the  case  of  finite  friction.  B.  Any 
set  of  loads  that  will  not  cause  collapse  when  all  the  coefficients  of  friction 
are  zero  will  not  cause  collapse  with  any  values  of  the  coefficients. 

Occasionally  these  theorems  enable  the  limit  load  to  be  computed 
precisely  for  finite  non-zero  friction.  The  well  Vnown  two-diinenslonal  footing 
nroblem  (which  will  be  diseussed  later)  for  a  von  Hiaes  material  providea  such 
an  example.  Two  solutions  are  available  for  upoer  bound  computations.  One 
by  Prandtl  (12)  contains  a  rigid  region  which  acts  as  an  extension  of  the 
punch  and  the  contact  area.  The  other  by  Hill  (13)  assumes  zero  friction,  and 
appreciable  slip  does  taVe  place.  Both  solutions  give  the  same  answer  for  the 
avenge  pressure  (2‘'’‘rr)  where  k  is  the  ’4eld  stress  in  simple  shear. 
This  value  also  has  been  shown  to  be  a  lower  bound.  Therefore,  the  limit 
nressure  is  k  for  all  possible  values  of  the  coefficient  of  friction. 

It  should  be  noted  that  while  the  limit  theorems  are  valid  for  an 
elastic-oerfectly  plastic  material,  the  elastic  strains  are  not  considered 
when  establishing  a  kinecatically  admissible  velocity  field.  Therefore,  the 
material  can  be  considered  to  be  rigid-perfectly  plastlci^  for  the  purpose 
of  determining  bounds  on  the  collapse  load}  for  a  rigid-perfeotly  plastic 
material,  as  well  for  an  elastic-perfectly  plastic  material,  three  groups  of 
differential  equations  describe  the  behavior  of  the  material,  namely  the 
stress  equilibrium  equations,  the  stress-strain  relations,  and  the  strain- 
disolacenBnt  relationships.  In  general,  an  infinity  of  stress  states  will 
satisfy  the  stress  boundary  conditions,  the  equilibrium  equations,  and  the 
yield  criterion}  and  an  infinite  number  of  displacement  fields  compatible 
with  a  continuous  distortion  can  be  found,  independent  of  the  stresses,  which 
satisfy  the  disnlacement  boundary  conditions.  The  stress-strain  relations 


A  rigid-perfeotly  plastic  material  is  one  which  undergoes  no  strain 
OTtil  the  n  Xa  j.c  state  is  reached,  in  which  deformation  occurs  at 
a  constant  state  of  stress. 
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are  necessary  ;  delorrmne  sta'cs  of  stress  and  displacement  that  correspond. 

In  elasticity  theory,  a  unique  solution  of  stresses  and  diaplaceroents  exists, 
however,  there  is  not  unioue  solution  for  a  rigid-plastic  material.  The 
solutions  norrrally  presented  consist  of  a  defoncation  mode  valid  for  the 
entire  body  considered  and  a  corresponding  stress  state  for  the  region  con¬ 
taining  the  plaaticaily  defornving  material.  T^e  oompletsneaa  of  such  a  solu¬ 
tion  depends  on  whether  it  is  possible  to  find  a  stress  solution,  eatisfying 
the  boundiry  conditions  and  equilibrium  equations  for  the  region  assorted  to 
be  rigid,  -^at  does  not  violate  the  yield  ccxidltlona.  These  oartial  solutions 
have  been  referred  to  as  "incor^lete”  solutions  and  the  ones  for  which  tte 
stress  solution  has  been  extended  to  the  shol©  body  concerned  are  called 
’'oomolete"  solutions  (Hi). 

The  velocity  fields  associated  with  incomplete  solutions  are  used 
to  establish  kineratically  admissible  fields,  and  hence  they  fumi^ 
upper  bound  values.  The  complete  solutions  for  stress  and  deformations,  while 
they  are  indeed  acceptable  solutions  for  a  rigid  plastic  material,  cannot 
be  repaided  as  the  limit  to  the  solution  for  a  real  material  whose  rigidity 
constants  have  become  infinite.  Nevertheless,  the  collapse  loads  associated 
with  these  solutions  art  the  actual  collapse  loads  for  an  elastic  iraterial 
with  infinite  moduli  U-t^ld-plastic)  or  for  that  matter  finite  moduli  (elastic- 
plastic).  The  relationship  uf  the  deformation  modes  obtained  in  conplete 
solutions  to  those  for  an  elastic-plastic  niaterial  (lli)  is  unknown, 

A,li  Previous  Studies 

Several  studies  have  been  nade  which  illustrate  well  the  application 
of  the  limit  ■yieorerns  to  the  determination  of  bounds  on  the  collapse  load  of 
a  footing.  In  these  studies,  the  weight  of  the  material  has  been  neglected 
in  the  stress  equilibrium  equations.  Investigations  regarding  the  influence 
of  the  weight  of  the  mass  have  not  oassed  beyond  the  stage  of  establishing 
the  differential  ecruatlons  (1).  The  weight  of  the  material  comnlioates  the 
situation  considerably.  At  given  values  of  cohesion  and  angle  of  internal 
friction,  the  material  weight  increases  the  critical  load  and  changes  the 
shape  of  the  surfaces  of  sliding  within  the  plastic  regions  of  the  material. 

The  problem  of  comnutlng  the  critical  lead  on  the  assumption  that  the  unit 
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waight  is  greater  than  gero  has  been  solved  only  by  aoproximata  methods  (1). 

The  above  mentioned  atudLes  are  baaed  on  the  aolutlons  prooosed  by 
Prandtl  and  Hill  or  generallgations  thereof.  Prandtl  (12)  proposed  a  stress 
field  at  inoipient  flow  for  the  plane  strain  problem  of  a  flat  rtgid  punch  on 
a  serai-infinite  mass  of  rigid-perfectly  plastic  material  which  obeys  Couloitb's 
yield  functioni  The  stress  field  was  given  only  for  the  plastic  oortlon  of 
the  mass  in  the  immediate  vicinity  of  the  ounehi  No  stress  solution  was 
oresented  for  the  rigid  region  away  from  the  ouneh.  The  velocity  field 
associated  with  the  stress  field  implies  that  the  punch  was  rough,  that 
no  relative  motion  occurs  between  the  punch  and  the  contact  area.  The  bearing 
pressure  corresponding  to  the  stress  field  is 

p  D  G  cot  0  0  ^  tan^  ( ^-  +  — ^— )  “  ^  ^  »  (A. 2) 

where  c  la  the  cohesion  and  0  is  the  angle  of  Internal  friction  in  Coulomb' 
law. 

Hill  (13)  later  also  oroposed  a  stress  solution  for  the  same 
oroblem  for  the  soeoial  case  where  0  •  0  (Tresca  or  von  Maos  function). 

This  stress  field  resulted  in  the  same  bearing  pressure  as  for  the  Prandtl 
solution]  this  nreasure  is 

p  =  (2  -hTT)  k  ,  (A.3) 

where  k  is  the  yield  stress  in  simple  ^ear.  Note  that  Zq.A.2  reduces 
by  use  of  L'Hospitals  rule  to  Eq  A.3  when  0  »  0  ,  where  c  corresponds 
to  k  .  The  velocity  field  associated  with  Hill's  solution  laplies  that  the 
ounch  is  smooth. 

Both  of  these  solutions  are  "incomplete"  solutions  as  defined 
earlier,  and  their  presentation  occurred  orior  to  the  development  of  the 
uoper  and  lower  bound  theorems  of  ll.-uLt  design.  However,  after  publication 
of  the  limit  theorems,  it  was  shown  that  the  velocity  fields  associated  with 
Hill's  and  prandtl's  solutions  were  kinematically  adnissible  fields  (l5). 
Therefore,  the  pressure  value  in  Eq3.A.2  and  A.3  are  upper  bounds  for 
Coulomb  and  Tresca  materials  respectively.  In  addition.  Shield  and  Itucker 
(15)  used  limit  analysis  to  show  that  a  lower  bound  on  the  collapse  pressure 
for  a  Tresca  material  is  5  k  . 
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More  recently,  Shield  (16)  has  extended  the  Prandtl  stress  solution 
into  the  rigid  region  and  succeeded  in  establishing  a  statically  admissible 
stress  field  throughout  the  material  for  values  of  0  <  75“  •  Therefore, 
the  Prandtl  value, Eq,  A. 2, is  both  an  uoper  and  lower  bound  on  the  collapse 
pressure  for  values  of  0  <  0  •<  75®  ,  and  hence  is  the  exact  value  of 
average  collaose  nrossuro. 

Shield  and  Drucker  (15)  extended  tho  application  of  limit  analysis 
to  the  three-dimensional  punch  oroblem  by  considerinB  the  problem  of  a 
rectangular  flat  nunch  on  a  Tresca  material.  It  was  shown  that  a  lower  bound 
for  any  rectangular  punch  is  again  5  k  while  the  upper  bound  for  a  siwoth 
Dunoh  lies  between  5*71  k  for  a  square  and  (2  +7T)  k  for  a  vary  long 
rectangle . 

Shield  (5)  subsequently -constructed  a  yield  surface  for  three- 
dimensional  stress  fields  baaed  on  Coulomb's  yield  function  and  used  the 
lower  bound  theorem  of  limit  analysis  to  determine  a  lower  bound  for  the  bearing 
capacity  of  a  smooth  or  rough  rectangular  footing  on  a  soil.  The  stress  field 
was  an  adaptation  of  the  statically  admissible  stress  field  previously  used 
by  Shield  and  Drucker  (15)  for  a  Tresca  material.  The  lower  bound  obtained 
for  the  limit  oressure  was 

p  ”  c  tan^  (-^  -|-)  -t-  sin  0  sin^  0  4-(l  4.  sin  0)  <lj'4-sln^  0) 

+  2  c  tan  (“^  (A.U) 

The  limit  theorems  have  also  been  used  in  problems  related  to  the 
classical  ounch  oroblem,  that. is,  in  the  problem  of  the  plastic  Indentation 
of  a’ layer  by  a  footlngo  i  comolete  plane  strain  solution  to  the  problem 
of  elastic  flow  in  a  thin  sheet  of  perfectly  plastic  material  laid  on  a 
rough  base  and  compressed  by  a  smooth  flat  footing  has  already  been  obtained 
(17).  In  addition,  through  e.mnlo/r.ent  of  the  limit  theorems,  Scield  (17) 
determined  unoer  and  lower  bounds  on  the  average  indentation  pressure  for 
the  same  oroblem  extended  to  three  dimensions,  for  a  square  and  circular  footing. 
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Another  promising  use  of  the  limit  design  theorems  in  this  area, 
is  In  problems  with  inhomogeneous  media  (6).  They  are  apt  to  be  of  greatest 
use  in  such  oases  because  of  the  enormous  difficulty  in  obtaining  exact 
solutions • 

A, 5  Illustrative  Problem 

To  illustrate  the  application  of  the  limit  design  theorems  to  soil 
mechanics,  the  upoer  and  lower  bound  determination  of  the  critical  height 
of  a  vertical  unsuoported  bank  of  soil  is  presented  here.  A  statically 
admissible  discontinuous  equilibrium  solution  (8)  is  shown  in  Fig.  A.l, 
where  w  denotes  unit  weight.  As  mentioned  earlier,  a  statically  admissible 
stress  field  must  satisfy  the  equilibrium  equations,  the  yield  criterion, 
and  the  stress  boundary  conditions,  In  addition,  since  the  proposed  stress 
field  is  discontinuous,  the  direct  stress  normal  to  the  discontinuity  line 
and  the  shear  stress  narallel  to  the  discontinuity  line,  must  be  continuous 
across  it.  It  is  easily  seen  from  Fig,  A.l  that  these  continuity  requirements 
and  also  the  stress  boundary  conditions  (naraly,  stress  free  boundaries) 
are  satisfied.  It  remains  to  be  shown  then,  that  the  stress  field  satisfies 
equilibrium  and  the  yield  criterion  in  each  of  the  three  regions  shown  in 
Fig.  A.l. 

The  equilibrium  equations  are 


2£L  i.  ■  alSL  a  0 

S  X  3  y 

ilsu+_lfz  .-V  . 

^y 


(A.5) 


The  first  equation  is  identically  satisfied  in  each  region.  The  first  term 
in  the  second  equation  is  zero,  and  the  second  term  is  equal-  to  -w  for  all 
three  regions.  Therefore,  both  equations  are  satisfied  throughout  the  entire 
body. 


*»HOUB  BSS6ABCH  COUNDATION  OP  I  i.  1.  I  N  O  I  S  INSTITUTE  OF  TECHNOLOGY 


ARF  Project  No.  8193-7 
Interim  Reoort 


Coulomb's  yield  criterion  my  be  expressed  as  follows: 


C  so  -  crtan  0  ,  (A. 6} 

where  c  is  the  cohesion)  0  is  the  angle  of  internal  friction)  U  and 
(taken  aa  positive  for  tension)  are  the  shearing  and  normal  stress  respectively 
on  a  failure  surface.  Alternatively,  the  failure  criterion  may  be  espreaesd 
by 


R  S  e  cos  0  - 


Oxt-  cn 


^  sin  0 


(A. 7) 


where  R  is  the  naxJ.mum  shearing  stress  at  a  point.  The  maximum  shearing 
stress  in  Region  I  is  w  H/2  and  occurs  at  the  lower  ground  level.  In  Region 
II,  the  maximum  shearing  stress  is  equal  at  every  point  and  is  also  w  H/2, 

The  plane  maximum  shearing  stress  is  zero  everywhere  in  Region  m.  There¬ 
fore,  the  yield  criterion  is  satisfied  in  Region  I  If 


w  FI 


0  cos  0 


-w  H 

~-T- 


sln  0 


or 


_2. 

w 


c  cos  0 
(1  -  sin  0) 


(A. 8) 


In  Region  II,  the  following  inequality  must  hold  for  all  values  of  y  , 

w  H  H  -  Zvy  -t-  w  H  _j  M 

— g-  sc  cos  0  -  ■  ■  *  ^  —  ■  sin  0 


or 


u  H 

(  1  sin  0)  -  w  y  sin  0  ^  c  cos  0  . 

Note  that  if  this  inequality  is  valid  for  y  =  H  ,  then  it  is  valid  for  all 
values  of  y  in  Region  II.  Therefore,  setting  y  =  H  the  above  inequality 
may  be  rewritten  as 


„  2  c  cos  0 

~  w  (1  -  sin  0)  ’ 


(A. 9) 


ARHOUft  RR5SARCH  FOUNDATION  OF  I  NOIS  INSTITUTE  OF  TECHNOLOGY 


ARF  Project  No.  8193*7 
Interim  Reoort 


which  IS  identical  to  Eo. A.R.  This  may  be  written  as 


tan  .  (A. 10) 

In  Ref^ion  ITI,  the  condition  to  be  met  is  that 

Oasc  cos  0  -  ^  y  sin  0 


or 


(H  -  y)  cot  0  , 

which  is  satisfied  for  all  values  of  y  in  Region  HI;  hence,  this  inequality 
nrovides  no  real  restriction  on  H  < 

Therefore,  if  the  condition  of  Eq.A.lO  is  satisfied,  the  stress 
field  in  Fig,A«l  is  statically  admissible  and  a  lower  bound  on  the  critical 
height  is  given  by  Eq.A.lO. 

A  velocity  field  with  a  rigid-body  slide  discontinuity  is  shown 
in  Fig.  A, 2.  We  shall  now  nroceed  to  orove  that  it  is  kinematically  admissible 
One  of  the  requirements  for  a  kinematically  admissible  field  was  that  the 
velocity  component  be  zero  on  those  portions  of  the  surface  where  the 
corresponding  component  of  the  surface  traction  is  not  prescribed.  However, 
in  this  problem  the  surface  tractions  are  prescribed  on  the  entire  surface, 
and  therefore  this  requirement  is  automatically  satisfied. 

Furthermore,  in  previous  discussion,  it  was  stated  that  for  a 
discontinuity  in  tangential  velocity,  such  as  is  proposed  in  Fig.  A. 2,  the 
angle  between  the  resultant  relative  velocity  vector  and  the  tH scontinuity 
surface  is  equal  to  0  ,  the  angle  of  internal  friction.  Therefore,  the 
normal  and  tangential  velocities  must  be  related  by 

3v'  -  3u'  tan  0  ,  (A. 11) 

as  is  shown  in  Fig.  A. 2,  where  '5v'  and  Su'  denote  the  normal  and 
tangential  velocities  respectively. 
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The  only  remaining  requirement  for  a  kinematically  admissible 
velocity  field  is  that  the  rate  at  which  the  actual  surface  tractions  and 
bocfy  forces  do  work  baaed  on  the  proposed  velocity  field,  must  equal  or  exceed 
the  rate  of  dissination  of  internal  energy  comnuted  from  the  strain  rates 
treated  as  purely  plastic  strain  rates. 

The  surface  tractions  do  no  work  since  they  are  prescribed  to  be 
zero.  The  rate  at  which  the  body  forces  (weight)  do  work  is 


w  ^3u'  coa^  -  3v'  sin^^  , 

From  Eq, A.lljthis  can  be  written  in  the  form 

w  t&njS  co^  -  tan  0  sir^  ,  (A.  12) 

Dissipation  of  internal  energy  takes  olace  only  at  the  surface 
of  sliding  since  the  sliding  body  is  rigid.  The  rate  of  dissipation  of  energy 
per  unit  area  of  surface  is  given  by  the  expression 

D  -  TT^u'  -hCrav'  ,  (A, 13) 

With  Eq,A,6  and  A, 11  this  equation  reduces  to 

D  »  c  du'  .  (A.ll;) 

The  total  rate  of  dissination  is  obtained  by  multiplying  D  by  the  length 
of  the  line  of  discontinuity,  H/cos^  ,  therefore 


c  ♦u’  E/co3yS  (A.l5) 

is  the  total  rate.  The  requirement  for  a  kinematically  admissible  field 
is  that  the  expression  in  E^,A.12  must  be  greater  or  equal  to  the  expression 
in  Eo, A,l5,  hence  ^ 

— V  H  tan^d  ^cos^  -  tan  0  sin/^  Ec/oos^ 
or 

— ~  w  H  sin/5  Tcos  0  cos/^  -  sin  0  sin^  -c  cos  0^0  . 
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This  reduces  to 


H  -s**  2  c  cos  0 

w  sin  /S  cos  (0+/S) 

Kdnimizlng  the  right-hand  side  gives 

Therefore, 

JJ  ^  _ ; _ 2  e  cos  0 _ 

w  sin  ( — ^  *  “2  ' )  cos 


(A.  16) 


(A. 17)  ' 

H  K.  (-^  -f  -|-) 

orovides  an  uoner  bound  on  the  critical  height,  .  Note,  from  Eq.A.lO 
and  A, 17,  there  is  a  factor  of  two  betHu-en  the  upoor.and  lower  bound, 

2^  tan  (.Jll  +  ^)  Li  tan  (-f  +  4-  )  •  ‘(A-lS) 

a 

The  upper  bound  can  be  improved  by  considering  a  rotational 
discontinuity  (logarithmic  spiral)  Instead  of  the  translational  type, 
however  the  translational  type  equal]y  wall  illustrates  the  method, 

A. 6  Dynamie  Problem 

Although  the  limit  theorems  provide  a  very  useful  tool  for  determining 
bounds  on  the  static  collapse  load,  no  such  general  methods  are  available  for 
the  case  of  dynamic  loading.  A  numerical  procedure  can  of  course  be  employed, 
however  the  addition  of  inertia  terms  in  the  governing  equations  will  certainly 
add  to  the  comnlexlty  of  the  oroblem  and  may  necessitate  more  simplifying 
assunntions  than  introduced  in  the  static  oroblem.  The  plane  plastic  strain 
oroblem  with  inertial  efft^cts  is  discussed  in  detail  in  the  following 
appendix. 


ASMOUB  BESEARCH  EOUNOATION  Of  ILLINOIS  INSTITUTE  OE  TECHNOLOGY 


ARF  Project  No.  8193-7 
Interim  Report 


A, 7  References 

1.  K.  Terzaghi,  Theoretical  Soil  Mechanics,  John  Ifiley  and  Sons,  19li3t 

2.  r>.  C.  Drucker,  R.  E.  Qihaon,  and  D.  J.  Henkel, "Soil  Mechanics  and 
Wurk-Hardening  Theories  of  Plasticity”,  Proo*  ASCE,  7ol.  61, 

Paoer  798,  Sopt,  1955. 

3.  W.  Prager  and  P.  G.  Hodge,  Theory  of  Perfectly  Plastic  Solids, 

John  Wiley  and  Sons,  1951. 

Ii.  D.  C.  Drucker,  ’’Some  Inrolications  of  Work-Hardening  and 
Ideal  Plasticity",  Q.  Appl,  Math..  7  (1950),  mi-ll8. 

5.  R.  T.  Shield,  "On  Coulomb's  Law  of  Failure  in  Soils", 

J,  Heeh.  Phya,  Solids,  h  (1955),  10-16. 

6.  D.  C.  Drucker,  "Limit  Analysis  of  Two-  and  Three- DLmenaional 

Soil  Mechanics  Problems",  J.  Heoh.  Phys.  Solids,  1  (1953),  217-226. 

7.  R.  T.  Shield,  "Mixed  Boundary  Value  Problems  In  Soil  Mechanics", 

J.  Hech.  Phys.  Solids,  11  (1953),  61-75. 

8.  D.  C.  Drucker  and  W.  Prager,  "Soil  Mechanics  and  Plastic  Analysis 
or  Limit  Design",  Q.  Appl.  Math.,  10  (1952),  157-165. 

9.  D.  C,  Drucker,  W.  Prager,  and  H.  J.  Greenberg,  "Extended  Lindt 
Design  Theorems  for  Continuous  Madia",  Q,  Appl.  Hath. ,  9  (1952), 
381-389. 

10.  R.  T.  Shield,  "Stress  and  Velocity  Fields  in  Soil  Mechanics", 

J,  Math.  Phys.  33.  2  (195Ii),  llili-l56. 

11.  D.  C.  Drooker,  "Coulomb  Friction,  Plasticity,  and  Limit  Loads", 

J.  Appl.  Mech.,  21  (195Ii),  71-7li. 

»♦  .» 

12.  L.  Prandtl,  "Ober  die  Harte  olastischer  Koroer",  Goettinger  Nachr. , 
Math.-Phys.  Kl.,  (1920),  7L-85. 

13..  R.  Hill,  "The  Plastic  Yielding  of  Notched  Bars  Under  Tension", 

Q.  J.  Hech.  Appl,  Hath..  2  (19li9),  bO-52. 

lU.  J.  F.  W.  Bishop,  "On  the  Complete  Solution  to  Problems  of  Dafona- 
tion  of  a  Plastic  Rigid  Material",  J.  Hech.  Phys.  Solids,  2  (1953), 

h3-53. 

15.  R.  T.  Shield,  and  D.  C.  Drucker,  "The  Apnlication  of  Limit  Analysis 
to  Punch-Indentation  Problems",  J,  Appl.  Mech.,  20  (1953),  h$3~h60. 

16,  R.  T.  Siield,  "Plastic  Potential  Theory  and  Prandtl  Bearing  Capacity 
Solution",  J.  Appl.  Mech.,  21  (195U),  193“191a. 


17. 

A  R  n  o  u  R 


R.  T.  Shield,  "The  Plastic  Indentation  of  a  Layer  by  a  Flat  Punch", 

R  *1  LLi  HO°%  INSTITUTE  OP  TECHNOLOGY 


ARF  Project  No.  8193-7 
Interim  Report 


1 

h 

[ 

Region  Ii  0^  -  *7^  ■ 

(Ty-  -Wy 

Region  ITltC  •  (r< 

^  JT 

•My  -  h: 

Region  III  ‘ 

T  “  ( 
*y  ' 

iTy-  -v<7 

e 

■ 

o 

Pig.  A.l  AW  ZQTOTBRIUH  SOLUTTOH 


Pig,  A,2  A  VELOCITT  SOLPTTON 


ARHOUB  RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOOT 


ART  Project  No.  8193-7 
Interim  Report 


Aooendlx  B 


THE  METHOD  OF  CHARACTERISTICS 
APPLIED  TO  PROBLEMS  OF  PLAKE  PLASTIC  STRAIN 
WITH  INERHAL  EFFECTS 


by  Phllin  0.  Hodge,  Jr, 
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Appendix  B 

BIS  HSTHOD  OF  CHARACTERISTICS 
APPLIEP  TO  PROBLEMS  OF  PIANE  PIASnC  STRAIN 
WITH  INERTIAL  EFFECTS 

by  Philip  0.  Hodge,  Jr. 


B.l  Introduction 

If  a  slowly  increasing  load  Is  applied  to  a  struobure  made  of  a 
rlftld-perfootly  plastic  naterial,  tho  striaaturo  will  romain  rigid  so  long 
as  the  load  remains  less  than  a  certain  critical  load  variously  l<nown  as 
the  "yield-point  load",  "oollapBe  load",  or  "limit  load".  If  tho  yield- 
point  load  la  reached  and  maintained  for  a  finite  length  of  tine,  the 
struoturo  will  deform  in  a  ouasi-statle  flow,  which  will  continue  indefinitely 
or  until  the  geometry  of  the  struoturo  is  auffielently  changed  so  that  the 
yield-point  load  of  the  deformed  struoture  is  different  from  ttiat  of  the  original 
structure.  However,  if  the  load  ie  irwreased  abo*/e  the  yield-point  load, 
there  will  bo  no  possible  equilibrium  configuration  of  etresses  and  all  or 
part  of  the  structure  will  engage  in  accelerated  plastic  flow.  In  conventional 
applications,  such  accelerated  plastic  flow  would  be  disastrous  to  the  struoture, 
and  its  theoretical  details  are  generally  of  little  practical  interest.  How¬ 
ever,  If  the  time  duration  of  the  overload  is  sufficiently  small,  then  the 
inertial  resistance  may  be  great  enough  to  limit  the  deformations  to  structurally 
reasonable  values,  Tho  relation  of  this  type  of  loading  to  blasts  from 
nuclear  or  H  S  detonations  Is  obvious. 

Previous  work  on  overloads  of  this  type  has  been  done  on  beam  and 
frame  type  structures,  (1-  circular  plates  (16  -  19)  and  circular 

cylindrical  shells  (20  -  26)  .  Tte  present  appendix  is  concerned  with 
problems  in  which  the  plastic  structure  is  in  a  state  cf  plane  strain.  A  typical 
problem  in  this  category  is  that  of  a  semi-inflnita  plastic  laass  pressed  upon 
by  a  rigid  punch.  (The  term  "structure"  is  perhaps  used  loosely  here,  re¬ 
ferring  as  it  does  to  the  semi-infinite  plastic  mass.)  Such  problems  have 


^  Numbers  in  parentheses  refer  to  the  references  listed  at  the  end  of  the 
appendix. 
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been  previously  trented  from  the  quasi-static  viewoolnt  (see,  for  ejsmple, 

27  or  28  for  text-book  treatments)  but  this  is  believed  to  be  the  first 
attenpt  to  inclucte  inertial  effects  in  a  nrobleni  of  olano  strain. 

In  the  next  section,  we  shall  first  state  the  basic  differential 
equations  of  the  oroblem  and  simplify  them  somewhat  by  a  change  of  variables, 
tfe  shall  then  find  the  cteraoteri sties  of  this  system  and  the  relations  that 
must  hold  along  those  characteristics.  This  discussion  will  be  carried  out 
for  doformationo  of  arbitrary  magnitude.  Section  B.3  will  begin  with  siraoli- 
flcation  of  the  eouitions  for  the  case  of  :;mall  deformations  such  that  the 
Eulerian  and  Laprangian  coordinates  need  not  be  distinguished  from  each 
other.  Iho  differential  equations  wil]  then  be  replaced  by  corresponding 
differencG  equations  in  prenaratton  for  n  numerical  solution.  The  following 
section  will  bo  onnoornod  v/ith  the  oxamole  mentioned  above  in  which  a  rigid 
ounch  is  nressed  against  a  sert. -infinite  plastic  mass.  A  possi.ble  numerical 
sohema  for  the  solution  of  this  example  will  be  presented,.  Finally,  the 
renort  will  conclude  with  some  suggestions  for  extension  of  the  results, 
both  to  other  oroMoms  of  nlane  strain  and  to  some  related  oroblenis  in  soil 
mechanics. 


D.  2  Basie  Equaf'^ns 

Under  the  usual  assumotions  of  nlane  strain,  three  stress 
comnonents,  ,  and  two  velocity  comnonents,  u.'  and 

y'  must  be  determined.  Primes  are  used  to  denote  physically  dimensioned 
quantities;  unorimed  symbols  vri.ll  nresently  be  introduced  for  corresponding 
dimensionless  quantities.  To  determine  these  quantities,  we  have  available 
two  equations  of  motion 
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the  yield  condition, 


r(a;',  t;,!  a  r,')*  +  4(T;,f  -  4k*  -  o 


and  three  flow-law  equations. 
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d  x' 
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(c) 

(B.l) 

(d) 
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®  'ty'  '  ax'  ■'  ST',, 

Here  is  an  unknown  scalar  function  that  represents  the  roagnllaide  of  the 

plasUo  strain-rate  tensor.  We  introduce  dlmenslonlaaa  quantities  defined  by 

X'=  Lx,  x/2-  ^/f'  v'o  yzK/p^  V 

(B.2) 

<r/»zkcr;,  (Ty'^Zkff;,  'r;^«Zk'r*y,  t'  -  l Jf~ f t 


where  l_  is  a  tynical  length  of  the  oroblem,  so  that  Eq  B.l  become 
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Next,  doflne  uJ  as  the  dimensionless  mean  normal  stress  and 
9  aa  the  angle  from  the  negative  y  axis  to  the  first  shear  direction 
(see  27) f  which  leads  to 


(f) 


0^  =  CS^*/2.K  ts  i_o  +  ^  &in  2© 
=  <S;'/2k  =  CO  t.in  2© 
=  'C-',,/2k  =  -  cOte  2© 


(B.i*) 


Equations  B.I4  identicall/  satisfy  Sq.  B.3>  At  the  same  time,  we  oonblne 
SqiB«3t  B*3d|  0|  and  f  so  as  to  eliminate  2.  •  The  resulting  set  of  four 
equations  for  the  UThnowns  to  »  ©  ,  vjl  ,  and  v  is 

CO^  ©,  coJb  2©  f  ©y  5»\0  2©  —  UU)|  “  vVy  ■=  U.^  (a) 


uOy  r  <b\n  2©  -  ©y  cots  2©  -  -  vv^  o  (b) 

(B.5) 


u. ,  © 


(c) 


u.^  c.C5%  2©  *r  u.y  t»\n  2©  V,  iiin  2©  -  Vy  coi2©  =0,  (j) 
where  subscriots  are  now  used  to  indicate  differentiation. 

The  explicit  time  dependence  of  the  problem  occurs  only  in  the 
right-hand  sides  of  Eq,B.5.  For  the  purpose  of  finding  the  characteristic 
curves  at  any  time  -t.  of  Eq«B.5  ,  we  may  tenporarily  regard  the  right-hand 
si<tes  as  known.  The  characteristic  curves  are  then  defined  as  those  curves 
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T  across  which  somo  derivatives  of  t*i  ,  ©  ,  u  ,  and/or  v  nay  be  dis¬ 
continuous]  although  the  functions  themselves  are  continuous.  If  to  is  known 
along  such  a  curve  P  ,  then  the  variation  dw  along  T  is  also  known. 
Therefore,  ttie  oartial  derivatives  of  to  on  V  must  satisfy 


tO^  dy  -►  to.^  dy  «  dtO  .  (a) 

Similar  relations  must  hold  for  the  other  variables « 

©,  dx  ©y  dy  «  ds  (b) 

dx  dy  »  d  tt  {c) 

V,  dx  Vy  (iy  =  J  V ,  (d) 


In  general,  the  eight  Eq.  9.5  and  B.6,  regarded  as  linear  algebraic 
equations  for  the  eight  derivatives  ,  Wy  ,  ®t  >  ®r  »  » 

,  Vy  will  yield  a  unique  solution.  However,  a  characteristic  curve 
is  defined  by  the  property  that  these  derivatives  are  not  all  unique,  but 
may  have  different  values  on  the  two  sides  of  the  curve.  The  condition  for 
•tola  lack  of  uniqueness  is  that  the  determinant  of  coefficients  vanishi 
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After  sone  raanieulation,  this  requirement  can  be  shown  to  be  eoui valent  to 
the  equation 

+  2  cos»  20  ^  ^  ,  , 

<dx,  \d)(i  1  (B.6) 


i: 


(^1  ' 

\dinl 


*12. 


t>m  2S 


Therefore,  Eq, B.8  represents  the  desired  differential  equation  for  the  characteris¬ 
tic  curves  of  Eq«B.$.  Solving  for  ,  we  see  that 


Ccit.  © 


(a) 

(B.9) 


dx 


^  =  -ton  B 


(b) 


are  each  double  characteristics  of  Eq*3*5i  we  denote  them  as  first  and  second 
characteristics,  resoectivaly.  In  view  of  the  geonetrio  interpretation  of 
Q  ,  the  characteristic  curves  are  everywhere  in  Uie  directions  of  the  lines 
of  orincinal  shearing  stress. 


The  orecedlng  results  are  ictentical  with  those  obtained  for  the 
corrasoonding  quasi-static  orobleni  (27  ,  28),  To  obtain  the  relations  which 
must  hold  along  the  characteristics,  we  refer  Eq,B.5  to  curvilinear  coordinates 
^  and  in  the  first  and  second  characteristic  directions,  respectively 
(Fig.  B.l).  Denoting  the  components  of  the  velocity  vector  by  ^  and 
in  the  ^  and  -n  directions,  respectively,  we  may  write 


c  a\n  ©  dx  "  co^  ©  dy  ; 
d'n  =  Cob  ©  dx  +  0  dv^; 


u-yubiriB  cob  © 

(B.IO) 

V  «  -/U.  cDb  0  b'Ci  © 


Substitution  of  Eq,  B.IO  into  S.5  together  with  some  elementary  manipulation 
leads  eventually  to  the  four  equations 

(a) 

(B.n) 

+  )j.  e,^  ^  o  (b) 
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(cO-0)^  =  >0©^)  -^©J 

(oi-*-  o')^  =  -*■  (d) 

The  terms  on  the  rlghVhand  slcfes  of  EqiB^le  and  d  are,  of  course,  the 
inertial  effects.  If  they  vanish,  Eq  B.ll  reduce  to  the  well-known  results 
for  qnasi-statlo  plasticity, 

B.3  Difference  Equations  for  Snail  Deformations 

If  the  total  deformations  are  sufficiently  small,  vs  need  no 
longer  distinguish  between  Eulerlan  and  Lagranglan  coordinates.  Ihus, 
regarding  the  left-hand  sides  of  Sq,  B.lb  and  c  as  written  in  Sulerian 
coordinates,  vs  may  use  the  siznple  time  derivate  of  the  middla  neuters 
rather  than  the  oxtreme  right-hand  sides  of  Eq,B.la  and  b.  Following  this 
argument  through,  wo  find  that  the  non-linear  velocity  terms  in  Eq,  B.llc  and 
d  have  disappeared  so  that  our  set  of  eouatlons  is  now 

^  j  ^  ®  O  (a) 


1ft  o 

(b) 

(oo  -  o)^  •= 

-  -^©^ 

(=) 

Cuo  -c  e)^  - 

(d) 

To  replace  Eq,B.12  by  difference  equations,  we  consider  first 
derivatives  with  respect  to  time.  At  any  instant  t  ,  we  may  reasonable 
assume  that  we  have  solved  the  problem  for  all  previous  times,  in 
particular,  at  a  time  -t -At.  •  If  we  excand  any  function  ^  at  -t-At.  in  a 
Taylor’ s  series  about  t  we  obtain 
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9  (t.  -  9t  ^  Ct)  £^'fc’‘  + 


•  «  «  » 


'IhersforB, 


p^(t)  =  1-^  +  k^t, 


(B.13) 


whore  k£^-t  represents  the  error  involved  In  truncating  the  series  and  the 
notation  cf  is  used  to  designate  , 

Replacing  the  thae  darlvates  In  Eq,B.l  by  Eq.Bil3*  ve  obtain 
the  set  of  equations 

(a 

(w-el  «  +  kAt  (0 


Cb.U) 


(0  -  ^)  4- 


k  ii't 


To  replace  derivatives  with  respect  to  ^  and  by  finite 
differences,  we  exoress  the  values  at  the  mesh  ooints  by  Taylor's  series 
about  the  raldoolnts  A  or  B  in  Fig,  B.2.  Thus 

<p.i 

J 

whence,  by  adding  or  subtracting. 


A^tj (a)  =  ^  kA.^\^ 


(B.15)  I 
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Similarly, 


9  Cb)  =  ^  (  ^ 


(c) 


(a) 


(B.16) 


Substitution  of  Eq, B.l?  into  Eq, H.lii  yields  the  four  difference  equations 

'^^•■4  ->*4  'N  *‘''=^5' 

-■^‘i  *t  ° 

(b) 

-  (  CO..  -  +■  -r  A  B. 

■“I  [(>^^-',4  ■*■ 

(co^j.,  -r  e^,.)  -  6tj')  -*-  k,  CA.  a.'T^  -».  kj  A'Tj^'^  «  (B.17) 

/  _  o 

In  the  anplicatlon  to  follov,  and  bj  will  generally  be  known 
at  points  to  the  left  and  below  and  It  will  be  necessary  to  oropagate  up  and 
to  the  rights  Therefore,  we  are  interested  in  solving  Eq,B.l6  for  and 

\)tj  in  terras  of  values  at  w-l  and  V.,  j-l  •  To  this  end,  we 
simply  replace  u  by  u-l  in  Eq,B.l6a  and  j  j"'  ^  Eq,B»l6b,  Doing 

this,  and  solving  the  resulting  linear  relations  explicitly  for  fjLi^  and 
O;.;  ,  we  obtain 


(B.18) 
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To  handle  the  terms  A  and  In  £q*  B.17  we  need  the 

geoaetrie  relatione  the  Cartesian  coordinates  Xi^  ,  ytj  •  By  use  of 
aporoortate  newer  series  exnansions.  It  Is  readily  verified  that 

°  i  A^ti(  a»>n  !>»n  ©tj) 

yu>.i  -V‘4  ^  urs»  +  k  A 

(B.20) 

(cos>0t,j„  ■*•  co^  ^n)  +  k  Ai^* 

v-i-  “  y^i  “  t  ®"'r‘ 

Solving  these  relations  for  x-'  i  'yij  >  A  ^^nd  Anr^- ,  we  obtain 

-  (tas,  'f  coa  ©ij)(yi..,^  -  yc,  ;*.)]  CB.21) 

-+-  (s.in  ■*■  S.j )(  ~  (^) 
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(B.22} 


>*4  =  -c  cos»  ©ij)(coi,  ©L,j*.  +•  ■*■ 

+  (am  ©ui.j  +  5>in  ©tj)(sin  ©{„j*i  ■'■  s»in  + 

0 

+  (t>tn  +  s.in  ©tOtcos,  9i,j„  +  cob  ©t/Xyc.,.}  -yt,i*.l] 

y^i  “  ^’•1  [(^"^  bmeij*)  + 

+  (  coi,  Qi„,j  +cob  ©4j  )(cob  at.,j,>  +  tob  ©tj)  -V- 

4  (cob  ■♦■  to6©nM*»'^  ®i.,;‘'  ■*"  ®‘j)(xc..4  -Xt,j*,^  . 

irtiore 

COb  -.©,.')+  cos,  (©t,j.,  ^ 

If  W8  subtract  Eq  B.6b  from  Eq  B.17a  and  use  Eq  B,21,  we  obtain  a  non-linear 
algebraic  equation  »rtiich  involves  only  S  at  the  point  I  ,  ^  i 
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(B.2U} 


It 


To  sommarlze,  ve  have  six  9qaatlon8(Eq,B.l8,  B.19>  B.2?,  B.23,  B.25  and  B.26) 
aTailabe  at  each  point  L,  j  to  determine  the  alx  variables  » 

Xij  I  yij  »  and  at  each  point.  We  shall  demonstrate 

a  possible  Iterative  mo'tiiod  for  solving  a  oortioular  problem  in  the  next  seetion. 

B.ii  ExMtple 

Ve  consider  as  an  exanple  a  rigid  ounoh  of  width  2.  L  pressed 
against  a  plastic  half-space  with  a  total  force  2  PK  •  The  qaasi-statia 
solution  (Fig.  B.3)  to  this  problem  is  veil  known  (see,  for  esampls  27).  If 
^e  punch  surface  is  oerfeotly  smooth,  then  there  exist  alternative  velocity 
fields  (27)^  To  avoid  possible  ambiguity,  wa\hall  assume  the  punch  surface 
to  be  rough. 

The  yield-point  load  is  P  »  ,  and  the  incipient  velocity 

field  oonsists  of  a  downward  motion  of  the  center  seotlon  ABA  of  arbitrary 
voloolty  4  ,  togethor  with  flow  along  ^e  eharaet»‘iatlee  of  taagnitade 
iy/yr  in  the  remaining  plaatio  regions. 

Ve  shall  assume  that  the  punch  displacement  x.  is  known  as  a 
function  of  tine  and  that 


t  (o')  »  4(0)  ■  o 

BCt)  O  (B.27) 

Z.  (t)  ts.  Be  <<  1  8er  all  t  >  T 

The  last  restriction  enables  us  to  use  the  Eulerian  form  of  the  equations  and 
to  neglect  the  motion  of  the  boundaries.  We  wish  to  find  the  solution  at  all 
times  t  and  in  particular  to  find  the  desired  force  P  as  a  function  of 
time. 

Ve  assume  that  the  solution  at  any  tins  b  has  the  general  character 
shown  in  Fig,  B.U.  A  central  curvilinear  wedge  AM  novas  rigidly  down  with 
velocity  a,  ;  in  the  region  ABC,  the  second  characteristics  all  pass  throng 
the  sLngular  point  Aj  the  characteristics  in  ACD  have  no  singular  poiat. 
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Along  the  wedge  boundary  A3,  the  normal  conDonent  of  velocity 
must  be  continuous,  hence 


B  ^  cos  &  on  A&  (B«27e) 

With  /Ji.  known,  Eo.  B.l^b  with  j  reolneed  bv  j-i  yields 

'Jy  »  -^  i  ©t^  ^B.27b) 

Since  0  must  vani^  at  B,  Eq.B.2  my  bo  SDolied  to  find  ^  at  each  point 
on  AB« 

On  tho  bound?ry  BCD,  the  norral  velocity  coraoonont  must  vanish. 


OaO  on  aCD 

It  then  follows  from  Eq,  DqIIh  that 


(B.28a) 


yo.  a  const. 


TlJ'JIL  on  QCD 


(B.28b) 


where  tho  constant  is  evaluated  from  continuity  of  /JU  at  B. 

On  AS,  no  loads  are  applied,  hence  it  follows  from  Eq  B.i*  that 


to  -  */z.  »  ® 


(B.29a,  b) 


Also,  since  motion  of  the  boundary  is  neglected 

y  eO  on  AD,  (B.29o) 

Finally,  if  we  prescribe  a  uniform  spacing  of  mesh  points  numbered  from  toe 
point  A, 


X),  S'  L  on  AD. 


(B.29d) 


It  is 

the  infinitesimal 
From  A ' '  to  A ' ' ' 
the  same  as  Bo  2?. 
counting  from  A' ’ 


convenient  to  regard  the  singular  point  A  e^janded  to 
curve  A'  A"  A'*-'  as  indicated  in  the  insert  to  Fig,  B.ii. 
the  field  is  not  singular  and  the  boundary  conditions  are 
However,  along  A'  A",  ©  on  the  -Lth  characteristic 
will  be  assigned  the  value 


©t  ~  (~^/4)  ^  A'^" 


CB.30a) 
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Of  course 


X  •  y  B  O  on 

so  that  0  and  E(j.Btl7a  redujes  to  to-0t*con8tant.  Hence, 


(B.30b,c) 


IaJ  a  -  ^  L  ^  om  A’  A»'* 

Assuming  uj  and  B  to  be  knoiisi  on  AB,  the  force  P  ie  given  by 


(BOOd) 


P  o  /  (-  lua  +  ■feair*  ©)  dx 


an  aquation  which  Is  easily  Integrated  numerically. 


(B.31) 


To  find  the  solution  at  any  time  t  ,  we  assume  It  known  at  time 
•te  -  At  <  Wo  guess  a  value  of  0  based  on  its  orevlous  values.  We  divide 
the  unknown  angle  between  the  chnractnrl sties  AB  and  AC  at  A  into  m 
equal  angles  ,  where  m  is  glvon  and  Is  guossod.  Similarly,  tho 

unknown  lenrth  AD  is  divided  into  a  known  number  of  Inten-nls  n  of  equal 
but  unknown  magnitude  Ax  •  Denote  the  point  C  by  0,  0  and  Increase 
along  CD  and  j  along  CA, 

Equations  3,27  and  B.28  show  that  /u  and  •?  along  AB  and  BC  are 
comoletely  suecifled  in  terns  of  the  guessed  values  of  ©  ,  Therefore, 
the  rl^t-hand  sides  of  Eq,D.l8  and  B.l?  are  all  known  and  and.  Oij  can 

be  found  at  points  next  to  the  boundaries  AB  and  BCD  •  Continuing  this 
process,  we  can  find  ^  and  0  throughout  the  field. 

We  next  apply  Eq  B,25  to  a  point  next  to  the  boundary  AD  ,  regarding 
everything  as  known  except  Qtj  .  Determining  a  corrected  value  for  Qcj  from 
the  resulting  non-linear  equation,  we  can  then  find 

Eq, B,22,  B.23  and  B.26,  respectively.  This  process  can  be  carried  out  success¬ 
ively  back  across  the  entire  field  and  eventually  furnishes  corrected  values 
of  e  in  the  boundaries  AB  and  ACD, 


This  new  solution  is  not  tested  against  three  criteria*  (1)  the 
value  of  ©  on  AB  should  not  ha\'e  changed  too  much  from  the  previous 
step,  (2)  ©  at  3  should  equal  'TT/U  ,  and  (3)  X  at  B  should  equal  -1. 
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Failure  of  either  or  both  of  the  last  two  criteria  calls  for  oorreoted  guesses 
for  ^  and  *  respectively,  and  failure  of  any  one  or  iwjre  of  the 
three  criteria  calls  for  an  iteration  of  the  entire  process^  using  the  new 
values  for  9  ,  A9  »  and  Ax  *  After  a  sufficient  nuinbep  of  iterations 
so  that  all  criteria  are  satisfied,  another  increment  of  time  la  added  and 
the  entire  process  is  repeated. 

B.5  Conclusions 

The  preceding  analysis  rooresonts  a  first  effort  In  i*at  is  be¬ 
lieved  to  be  a  new  field.  Some  imoortant  theoretical  questions  renaln  un¬ 
answered,  and  numerous  future  problems  are  suggested. 

Among  the  theoretical  questions  the  primary  ones  deal  wlti\ 
convergence,  uniqueness,  and  the  non-linear  terras  in  Eq,  B.llo  and  d. 

It  seems  Intuitively  plaueible  that  the  iteration  procedure  suggested  for  a 
particular  example  in  the  previous  section  will  converge.  However,  intuition 
is  no  substitute  for  nHthamatlcal  rigor,  and  it  would  be  most  desirable  to 
prove  this  eonvorgence  in  the  most  general  terms  possible. 

The  question  of  uniquenoss  is  always  a  thorny  one  for  a  perfectly 
plastic  material.  Indeed,  in  quasi-static  oroblems  only  the  yield-point  load 
is  unique,  it  frequently  being  possible  to  construct  alternative  stress  or 
velocity  fields.  To  the  author's  knowledge,  uniqueness  of  solutions  to  dynaolo 
plasticity  problems  has  received  no  attention  as  yet. 

Finally,  among  the  mathematical  ouestions  to  be  answered.  Is  the 
disturbing* one  concerning  the  non-linear  velocity  terms  In  the  equations  of 
motion.  It  will  be  recalled  that  the  characteristics  were  defined  as  those 
curves  across  which  some  of  the  derivatives  with  respect  to  might  be 
di^ontinuous.  Therefore,  It  is  difficult  to  accept  the  appearance  of  such 
derivatives  in  Eq,  B.llc.  The  whole  subject  of  the  nature  of  discontinuities 
permitted  across  the  characteristics  needs  to  be  more  fully  investigated  for 
equations  of  higher  than  second  order. 

Assuming  that  the  appearance  of  the  ^  terms  in  Eq, B.llc  is  satis¬ 
factorily  explained,  there  still  remains  the  question  of  how  to  evaluate 
the  equation.  Although  finite-difference  approximations  for  cross  derivative 
terms  are  easily  obtained,  they  apparently  give  a  lower  order  of  accuracy 
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for  the  same  number  of  ooints 


Tumlnp  next  to  other  related  problems,  an  entire  vista  opens  up. 
Indeed,  any  problem  for  which  the  quasi-statio  solution  is  icnown  would  seem 
to  bo  fair  pama  for  a  dynamic  analysis,  b.r,.  Indentation  of  a  plastic  mass  by 
a  ripid  wedpe,  Lookinp  still  further  ahead,  one  can  consider  rotationally 
syinmotric  nlastlc  oroblems  whose  quasi-static  solutions  have  been  obtained  by 
Shield  (29,  30). 

Still  another  direction,  and  one  of  considerable  oracttcal  iniaortance. 
Is  the  oxtension  to  tyic  failure  criteria  usually  used  in  soil  moohanies.  Hero 
the  yield  atross  k  ,  instead  of  boinp  constant,  la  a  function  of  oO  ,  Thus, 

Fq  n.lc  must  be  renlacod  by 

-  'O  (g.32) 

Obviously,  this  modification  tvill  affect  not  only  the  stress  equations,  but 
also  the  flow  law.  Hnupvcr,  althouph  the  rcsultihp  equations  will  certainly 
be  more  comnlicated  than  Eq,3.5,  they  should  be  of  the  same  order  and  contain 
the  same  derivatives.  It  therefore  annears  reasonable  to  hope  that  the  re- 
aultinp  oroblem,  nlthouph  more  comnlicated,  will  not  be  intrinsically  more 
difficult.  Since  even  for  the  nerfectly  plastic  material,  a  high-speed 
conroutlng  machine  must  be  resorted  to,  there  may  well  be  no  essentially  new 
oroblems  encountered  in  such  a  modified  analysis. 
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A-^oendix  C 


FOOTT^JDS  3!'\TECT2D  TO  ECCENTRIC  LOADS 


C.l  Introduction 

This  nppen<fl.x  connidera  footinEo  subject  te  eccentric  loads,  i.e., 
Qcmbinations  of  direct  vertical  loads  and  overturning  moments*  This  loading 
condition  la  ahovm  on  Pig,  C.l,  whore  tho  moment  is  Pe  ,  e  renroRenting  the 
eoul valent  eocontrlclty  of  tho  load.  Conalderatlon  of  the  noaslbility  of  a 
moment,  Po  ,  ineroases  the  neneraHty  of  the  solutions  oresentad  in  the 
Final  Henort^'^  on  tho  original  nropram.  It  ahould  ho  noted  that  since  the  ease 
where,  e  0  ,  reduces  to  the  ono-aided  failure  for  central  vertical  loads  as 
oonnldcred  oarllrr  ,  that  case  (o  »  0)  can  he  considered  as  q  seecial  class 
of  solutions  for  tho  rosults  nresentod  heroin, 

Tho  assumptions  made  in  this  annondlx  are  the  sans  as  those  made 
in  tho  oririnal  report  and  ronrooent  vet  a  further  extention  of  Andersen's 
method-'  . 

1.  The  failure  .'iurfaco  utidor  dyiiaMlo  loads  v.lll  be  tho  surface 
determined  by  nonlicatlor.  of  tho  initial  value  of  the  dynamic 
loads  as  a'  static  load. 

2.  The  resistance  offered  to  movement  is  rigid  plastic  in  form, 
i.e,,  settlement  and  soil  compressibility  are  not  considered. 

3.  The  maximum  clastic  resistance  eouals  the  static  resistance 
determined  analytically,  associated  with  the  failure  surfaces 
considered  in  item  1  above. 

h.  The  behavior  of  the  soil  is  defined  by  the  parameters  8", 

0  ,  and  c  ,  where  0  and  c  may  themselves  be  functions 
of  many  parameters  relating  to  the  soil  and  the  condition 
of  loading. 

First  the  static  oroblem  will  be  treated,  and  then  the  earlier  dynamic  anoroach 
will  be  modified  to  include  overturning  moments.  For  the  nurooses  of  this  aooendix 
pp.ly  the  reRUlting  ecuat'.or.s  will  be  nrecented  --  the  development  follows  that 
shown  in  earlier  reoorts  on  this  project. 

TT - : - 

-  McKee,  K.  E.,  Design  and  Ar^lvsis  oT  Feundntions  for  Protective  Structures, 

AFSWC,  TR  59-557 October,  ^ 

2/ 

-  Andersen,  P,,  Substructure  Analysis  and  Design,  the  Ronald  Press  Co., 

New  YorK,  New  York,  195^,  p.  t^l, 
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C.2  Static  Loads 


Based  on  Fip,  C.l,  the  static  solution  can  be  developed  for  a 
two-dimensional  footinp;  subjected  to  eccentric  loads,  following  Andersen's 
general  approach.  This  results  in  the  follovdng  two  equations!  I 


rr  ['■ECT  r  +.  D)  Da  -H  “  r | 
and 

a  r  3  *1 

LlJ«-  -  (2  -  ^  tan  jS)  -  -j-  (2  tan  -  1)+  (l-h—)^  (l-*-2  tan|) 

;  •  (0.2)'* 

[l+^tan  0j[^  H-(l  | 

Ihese  two  equations  are  coniparlable  to  Eq.  E.l  and  £.2  in  the  i 

Final  Report^/  •  The  first  aquation,  C.l,  is  identical  with  E.l,  while  the 
right-hand  side-of  C.2  is  identical  with  the  right-hand  side  of  E. 2.  Hence, 
the  only  Influence  of  the  eccentricity  is  on  the  left-hand  side  of  Eq.  C.2. 

This  allows  gross  simplification  malcing  use  of  the  earlier  work.  Certainly 
if  6  =  0,  the  eouatlons  are  the  same  and  Appendix  is  applicable. 

Furthermore,  when  e  «■  0  the  two  sets  of  eouations  are  identical  if  the 
quantity  (B  -  2e)  is  used  in  place  of  B  in  the  earlier  work.  With  this  , 

substitution  in  the  abscissa,  Fig,  E,2  through  E,26  are  apolioable.  As 

( 

before,  the  usefulness  of  the  plotted  solutions -is  limited  because  it  is 
not  practical  to  clot  for  all  values  of  q  . 

Equations  C.l  and  C.2,  of  course,  can  be  solved  num.erioally  for  ^ 

any  two  parameters.  For  the  applications  of  this  project,  the  paranaters 
that  would  ordinarily  be  determined  are  the  ultimate  load  capacity,  F  ,  ! 

®  I 

and  the  radius  of  the  failure  surface,  r  .  The  other  oarametersi  soil 

T7'  ""  "  •  "  '  •  - - -  ' 

McKee,  K.  E.,  Op.  Cit. 
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nrCDcrties,  0  and  c  j  footinp  dimensions,  B  and  D  ;  and  eccentricity, 

B  ,  would  have  to  bo  ooecified.  Trial  and  error  solutions  are,  in  peneral, 

necessary.  The  simolest  mothod  of  solution  rooidros  the  assumntion  of  a 

value  for  r  ,  which  allows  direct  solution  of  Eq,C.l  for  .  The  assumed 

r  and  the  commuted  P  can  then  be  substituted  into  Eq,  C.2  to  solve  for  B  . 

s 

This  Broeedm*0  is  reneated  as  often  as  required  until  the  desired  value  of  3 
is  determined,  and,  hence,  the  approorlate  values  of  and  r  are  establishod. 
Extraoolation  or  interpolation  can  he  used,  of  course,  to  find  values  when 
nearby  values  are  known. 

In  tho  oriplnal  work  for  e  o  0  ,  the  effect  of  q  was  con¬ 
sidered.  The  additional  factor  being  considered  here  in  the  eccentricity,  e, 
Figure  C.2  and  C.3  show  tho  Influence  of  e  on  and  r  resoectivaly, 

C.3  Dynamic  Loads 

Consider  tho  basic  dynamic  equation  asi 
10  4-  R  (9)  °  HCt)  ,  (C.3) 

where 


Q  ^  rotation  of  aoi 1 

..  ^  a 

Q  °  ^  acceleration 

dt'^ 

I  =  rotational  inertia 

R(Q)  =  resistance  as  a  function  of  0 

M(t)  =  time- depen dent  moments. 

For  the  faliure  natcern  anpreorinte  'o  fr.j  static  luadir.,:;,  ^hc  terms  in 
the  dynamic  ecuatton,  Eq,  C.3,  are; 


g(l  -f  ■ 


X 


0.6U357  r^-h  r^)-¥-  0.78510  r  -  D)^ 


(C.l 

ARMOUR  RtSEARCH  POUNDATION  Of  rLLINOIS  INSTITUTE  OF  TECHNOLOGT 


ARF  Project  No.  8l?3-7 
Interim  Heoort 


where 


R(©>= 


y© 

(r  +.0  -  -|-) 


(r  +•  D)^  tanvj 


D  (r  -  B)^ 
2 


(C.?) 


M(t)  n  P(t)  -  Pg 


0.106  YTC 

(r  4-0  -  4—) 


(r  -f  D) 


-h 


yp  (r  -  3)^ 
2(r  4-  0  '  -~) 


(C.6) 


y  o  unit  WQlKht  of  soil 

y-  tan"^  1^1  -  dJ 

g  °  gravitational  eonatant 

B  °  width  of  footing 

r  “  dlrfiension  defining  location  of  falluro  surface 

D  n  dooth  of  burial  of  footing 

9  '  rotation  of  soil  mass  about  ooint  C 


P(t)  ”  time- dependent  force  applied  to  footing 

P  "  static  caoacity  of  fo<.-ting. 
s 

The  above  developnwnt  uses  the  shear  surface  location,  r  ,  and 

the  load  caoacity,  P  ,  determined  by  the  static  analysis  nresented  in  the 

s 

previous  section.  It  should  be  emphasized,  however,  that  the  determination 
of  the  static  data  should  be  based  on  the  best  possible  estimate  for  the  soil 
orooertiesj  i.e.,  the  prooerties  should  be  sBlBct.eri  to  incorporate  the  influ¬ 
ence  of  the  variables. involved,  e.g,,  the  rate  of  load  application. 
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Appendix  D 

PROPSRTIES  OF  OTTAWA  SAim 


Host  of  the  pxporlmontal  work  of  the  original  and  ourront  programs 
has  been  conducted  using  Ottawa  saraS,  Oi»  of  the  puasling  aspaets  of  the 
original  research  was  the  dloeropaneloe  rested  between  the  properties  of  the 
Ottawa  sand  found  by  triasdal  teats  and  those  detamlned  by  calculations  based 
on  thn  footing  behavior^ •  This  appendix  contains  the  results  of  tvre  tests 
which  may  serve  to  clarify  this  situation. 

The  Ottawa  sand  ie  descrlbsd  by  the  grain  size  distribution  curve 
included  as  Fig,  D.l,  (The  original  sand  Is  still  used  in  the  sand  box, 
and  a  new  sand  is  used  in  the  glass  box.)  Figure  D.2  shows  the  results 
of  direct  shear  tests  on  this  sand,  when  the  density  ranged  from  110,0  to  117«5 
pef.  (These  densities  nro  glvon  as  a  range  althou^,  due  to  the  small  speol- 
men  size  and  other  factors,  they  are  individually  not  too  necurato.  In  parti¬ 
cular,  the  upper  bound  is  higher  than  that  previously  obtained  and  must  ba 
quostionod.)  The  envelopes  for  the  maximum  and  minimum  values  of  the  load 
deflection  curve  are  indicated. 

A  series  of  trlaxial  tests  wore  conducted  using  specimens  prepared 
in  the  harvard  Miniature  Compaction  liold  (1, 31-ln,  diameter,  2.82-ln,  height), 

A  number  of  alternate  approaches  for  obtaining  high  density  were  investigated. 
It  was  found  that  densities  of  approximately  10?  pcf  could  be  obtained  by 
compacting  the  sand  in  six  layers  using  a  weight  of  100  dropped  50  times 
per  layer  from  a  height  of  1  in.  For  each  test,  a  value  of  the  angle  of 
internal  fraction  0,  was  determined  by  passing  a  line  through  the  origin 
tangent  to  the  Hohr's  circle.  Table  D.l  and  Fig.  D.3  show  the  resulting  data. 
The  curve  and  points  on  Fig.  D.3  give  an  indication  of  the  scatter  In  the  data. 

The  experimental  difficulty  arises  here  in  the  determination  of  the 
density.  For  the  mnld  being  used,  a  l-gm  difference  in  weight  represents  1  pcf, 
In  addition,  the  volume  of  the  trlaxial  specimen  at  the  time  of  testing  may 
differ  from  the  volume  of  tne  mold.  With  respect  to  these  factors,  variations 

^ McKee,  K.F,,  Design  and  Analysis  of  Foundations  for  Protective  Ftraotures, 
AFSWC  TR  59-56,  October,  1939. 
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of  1  pcf  from  the  densities  reported  seen  reasonable.  In  spite  of  these  in¬ 
consistencies,  the  results  do  indicate  that  0  increases  8v4)8tantlally 
with  density. 

The  results  of  those  two  teat  series  ^*2  and  D,3)  establish 

that  the  angle  of  internal  friction,  0  ,  for  danse  dry  Ottawa  sand  can  exceed 
ho*  and  that  the  exact  \’alu0  is  a  function  of  the  density.  The  value  of  0 
found  by  inverse  application  of  Torzarhi'a  fonnula^  therefojre  may  represent 
an  accurate  valuo  for  the  sand  denel  ty  encountered. 


^  McKee,  K.  K.,  Op.  Ctt. 
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Table  .1 


»  B  M  O  U  B 


RESULTS  OF  TRI AXIAL  TESTE  OU  OTTAVtt  SlUD 


Density,  pcf  Anpclo  of  Internal  Friction, 

_ _  _  deg. _ _ 


102.0 

26,5 

105.2 

31.0 

105.8 

33.0 

107.1 

3li.O 

107.5 

33.5 

107,7 

35.0 

107.7 

33.5 

107,9 

33.0 

106,0 

37.0 

108.6 

32.5 

109.0 
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Appendix  E 

TWO-DIHENFIONAL  STUDIES 


E.l  Introduction 

Two-dtmonslonal  studloa  are  being  cairlod  out  In  a  glaas-alded  soil 
container!  This  container  was  doalgned  and  built  during  the  aarllor  raoaaroh 
program  and  lo  completely  doaoribed  In  earlier  roport^'*^.  The  objective  of 
this  phase  of  the  experimental  study  Is  to  obtain  qualitative  Information 
regarding  the  Influonoe  of  soil  typo  on  fbotlngo  oubjooted  to  both  statlo  and 
dynamic  loadn.  Bneaueo  of  the  nature  of  this  objootlva.  It  la  aasontlal  that 
this  phase  of  the  exporlmontal  program  bo  eonplotad  at  tho  earllost  poealble 
date,  alneo  It  la  nntlc^patnd  that  the  theoretical  approach  will  bo  guided 
to  a  major  extent  by  tho  obsorvod  resulta. 

The  range  of  polls  considered  for  this  phase  inolude  dense  Ottawa 
sand,  loose  Ottawa  sand,  dense  California  sand,  and  several  typos  of  oohesive 
sollg.  The  experlmentfi  ooneldered  Include  static  and  (^nam^o  loading#  For  the 
static  experiments,  the  procedures  are  followlnr  those  used  during  tho  original 
program.  For  the  dyoamlo  experiments,  tho  majcr  observations  are  being  made 
through  use  of  Faxtnx  photography  during  the  tost  and  through  bofore-ajxi-after 
still  photography.  These  tochnlquoo  Intended  primarily  for  qualitative  infonna- 
tlon  are  being  supplemented  by  quantitative  measurements  within  the  limits  of  the 
exporlmontal  procedures.  The  limitations  of  this  "quantitative"  Information 
should  be  clearly  understood  —  it  can  be  meaningful  or  It  may  be  satisfactory 
only  for  tho  purposes  of  comparison. 

E.2  Static  Tests 

The  glass  box  static  test-'  to  date  have  considered  prlrwrily  dense 
and  loose  Ottawa  sand}  one  trial  run  was  also  made  with  a  oohesive  material. 

These  experiments  were  conducted  with  normal  loading  (i.e. ,  through  a 
hydraulic  jack  with  the  total  loading  occurring  within  5  min.)  and  with 
mechanically  controlled  loading  rates  (I.e,,  the  downward  speed  of  the 
top  of  the  proving  ring  was  limited  to  a  specified  rate  of  0,000^3  in. 

McKee.  K.  E. ,  Op,  Git, 
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per  minute).  A  jnovie  camera  wag  automatically  turned  on  at  T?-  to  .’jO-mln. 
intervals  to  take  6  to  10  frames  of  the  latter  tosts.  Because  of  dlfflcultie-' 
in  reading  the  pages  from  the  movie  film,  only  United  data  vrero  cbtainrJ 
from  the  four  tspts.  Additional  experiments,  v^ere  the  time  lapse  photograph 
is  limited  to  dial  pictures,  will  be  required, 

The  results  of  the  glass  box  tosto  on  Otta^>a  sand  are  given  cn 
Fig.  E.l  and  E,?  for  donse  and  looso  aand,  roopoetlTOly.  Teats  FI  through  F7 
refer  to  normal  loading  rate  load  tests  using  the  hydraulic  jack.  Tests  FA  through 
FD  refer  to  long  term  experiments  at  controlled  loading  rates.  Jn  every  eano,. 
the  average  density  in  the  box  is  showi  after  the  symbol  for  the  test.  The 
United  data  from  the  exporlrmnts  with  controllod  rate  of  loading  showg  results 
subatantially  in  agraomont  with  those  for  normally  loaded  footings.  Ono  observa" 
tion  from  those  eontrolled-londing  rate  tests  is  of  particular  interest'-  the 
load-defloetion  curro  was  not  omoeth.  As  the  top  of  the  proving  ring  was  dis¬ 
placed,  the  load  incfogeed  with  esEentially  no  inernaso  in  displacement  until 
there  was  euddenly  a  proBo  motion  with  an  afisoelatod  decrease  in  load.  This 
cycle  was  repeated  throughout  the  duration  of  the  loading, 

Wth  the  controlled  loading  ratn,  the  formation  and  subsequent 
trano forma tion  of  the  soil  mass  wore  observable  In  proater  detail  than 
previously.  Figure  B,3  and  B,h  contain  a  selected  series  of  photographs 
for  such  tests  on  dense  and  loose  sand,  respectively.  The  behavior  of  footings 
on  dense  sand  is  similar  to  that  observed  in  the  original  program.  The 
footings  on  loose  sand  are  typified  by  the  photographs  included.  Large 
displacements  occur  by  local  shear  fallure^^  with  indication  of  a  shear 
surface  appearing  only  at  relatively  large  deflections  (see  Fig,  E.lie). 

E,3  Dynamic  Tests 

Eight  dynamic  tests  (six  on  dense  Ottawa  sand  and  two  on  loose 
Ottawa  aand)  have  been  conducted  in  the  glass  box.  All  of  these  experiments 
made  use  of  the  dynamic  loading  device  constructed  for  this  program.  Instru¬ 
mentation  was  limited  to  Fastax  photography  of  the  footing  response.  These 
experiments  were  conducted  in  an  attempt  to  increase  the  knowledge  available 
regarding  the  behavior  of  footings  subjected  to  dynamic  loads.  It  wis  hoped 

^  Terzaghi,  K.,  Theoretical  Soil  Neebanics,  Wiley,  19h3. 
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that  this  qualitative  Infoimatlon  would  aid  In  improving  the  analytloal 
taohnlquee,  while  dependenaa  was  placed  on  the  threp-dlmenei  onal  teste  for 
quantitative  data, 

PlRuroe  E.9  and  E.IO  show  a  nories  of  photographo  taken  fron  the 
faatax  film  for  danso  and  looso  Ottawa  oand,  rsEpeetlvely,  Four  of  the  beet 
fnetnx  fllme  (two  donee  sand  and  tvn  loono  aand)  wore  cut  ord  splleed 
to  fom  a  p^ntorlal  roeord  of  tt»  ojcpnrlmonta,  A  pri  nt  of  this  film  will  be 
made  available  on  request  to  AFFWC, 

A  single  plaoB  box  tost  was  conduotad  ufllng  eoheelve  material. 

The  donolty  In  the  box  Mae  121,0  pef.  The  loado  woro  applied  statically 
through  a  hydraulic  Jaek,  Tho  load-dlsplneoront  ourvs  Is  shown  on  E,S| 
and  photographs  appear  on  Pig,  B,6,  E.7,  and  E,fl,  The  strength  of  thlB  soil 
wna  sufficient  so  that  tho  glass  side  of  tho  box  ruptured  prior  to  formulation 
of  a  shear  surface  in  the  poll. 
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(c) 


Fig.  e.L  static  footing  failuhe  on  loose, sand 
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(e) 

Fig.  E.U  SIATIC  FOOTING  FAILURE  ON  LOOSE  BAND  (oont) 
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Fig.  E.5  FOOTING  BmVIOR  ON  COHESIVE  SOIL 
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Pig.  B.6  G0HE3TVS  SOIL  IN  QU3S  BOX  BKFQRB  TEST 
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Pig.  S.7  COHESIVE  apn.  IN  0U33  BOX  IHJRIMQ  TS3T 


Fig.  S.8  COHESIVE  SOIL  M'TER  GU53  IN  BOX  FAILED 
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Fig.  K.IO  DYMMIC  TEST  OF  FDOTES  OE  LOOSE  Omia  Rflm 
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THREE- DIMENSIONAL  STUDIES 


F.l  Introduction 

Three- Araenaional  studies  are  being  carried  out  in  the  sand  box 
used  during  the  earlier  research  on  this  project.  This  box  is  square 

in  olan  and  3- ft  deep.  The  objective  of  this  chase  of  the  experimental 
study  la  to  obtain  ouantltative  Infornation  on  the  behavior  of  footings 
subjected  to  static  and  dynamic  loads.  During  the  original  orojeot, 
static  loads  and  dropned  weights  were  aoplied  to  footings  placed  on  dense  dry 
Ottawa  sand. 


The  soils  to  be  considered  in  the  sand  box  tests  are  being  limited 
to  dense  Ottawa  and  dense  California  sand.  Within  the  scope  of  the  present 
orogram,  consideration  of  cohesive  materials  is  thought  not  to  be  practical, 
as  the  control  required  in  laboratory  exceriments,  the  mixing  and  placing 
of  cohesive  materials  offers  extreme  difficultyj  even  for  very  small  spacl- 
mens,  the  maintenance  of  uniformity  of  density  and  ingredients  requires 
rather  extreme  precautions.  Certainly  the  size  of  the  three- dimensional' 
container,  i.e,,  U  ft  x  U  ft  x  3  ft  ,  precludes  its  use  for  cohesive 
materials  without  elaborate  mixing  and  placing  apparatus.  If  in  the  future 
experiments  using  cohesive  materials  are  planned,  it  is  strongly  recommended 
that  consideration  be  given  to  using  the  smallest  soil  specimen  size  consistent 
with  the  required  results. 

The  three-dimensional  exDeriments  were  intended  primarily  to  yield 
quantitative  data.  For  this  reason,  emphasis  has  been  placed  on  selecting 
the  prooer  gages,  recording  s:,^tems,  etc.  Visual  observations  and  nhotographic 
records  have  been  nade,  but  only  for  general  coverage  of  the  experiments. 


F .2  Static  Tests 

As  cart  of  the  original  orogram,  a  series  of  footings  on  dense 

dry  Ottawa  sand  were  subjected  to  stati:  uads.  The  results  of  these  tests 

were  included  in  Aopendix  A  of  the  Final  Reporti'^  .  Additional  static  exoeri- 

1/  McKee,  K.  E.,  'Design  and  Analysis  of  Foundations  for  Protective  Structures, 
AFSWC  TR  p9-56TTctoFerr*I?p7^  ^ 
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mnnts  on  Oftaua  sand  were  not  renuired.  and  hence  no  static  experiments  have 
been  conducted  to  date, 

static  teats  will  be  conducted  with  California  sand  (shipped  from 
Naval  Civil  Engineering  Laboratory  for  use  in  Ws  study).  The  a1.m  of  thesu 
planned  experiments  will  be  to  obtain  suf fj cient  data  to  compare  the  reeulta  with 
those  for  Ottawa  sand.  It  la  anticipated  that  this  correlation  can  be  achieved 
by  use  of  the  failure  theories  for  footings.  If  this  proves  correct,  only  a 
limited  number  of  experiments  will  be  required. 

Three-dimensional  static  experiments  using  loose  sand,  either 
Ottawa  or  California,  could  be  conducted.  Placement  of  the  sand  in  a  loose 
state  would  require  a  major  effort  for  each  experiment.  The  advantage  of 
three-dimensional  experiments  using  loose  sand  vrLU  be  weighed  against  the 
complexity  of  conducting  such  experiments, 

F,3  Dynamic  Teats 

The  dynamic  loading  apparatus  is  described  in  Appendix  0,  Much 
of  the  developmental  experimentation  for  this  device  was  conducted  with 
footings  placed  on  the  surface  of  the  sand  box.  Although  meaningful  records 
were  obtained  during  many  of  these  preliminary  experiments,  it  was  decided 
to  make  use  only  of  records  obtained  using  the  final  apparatus.  These  records 
are  in  the  form  of  oscillograms.  The  data  are  a  force- time  record  from  the 
force  washer  and  a  displacement-time  record  from  the  LVDT  (Linear  Variable 
Differential  Transformer).  Visual  observations  are  made  before  and  after 
each  loading.  To  date,  only  limited  results  are  available  and  no  attempt 
has  been  made  to  analyze  the  records  obtained.  Examples  of  these  records 
are  given  in  Appendix  G.  Figures  F.l  and  F.2  show  examples  of  footing 
failure  under  dynamlo  loads.  All  of  the  dynamic  experiments  have  been 
conducted  using  dense  Ottawa  sand.  More  detailed  coverage  of  these  dynamic 
results  will  be  presented  in  the  final  report  on  the  current  project. 
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Appendix  0 

DYMAMIO  LQADHra  APPARATUS 


An  apparatus  suitable  for  applyinf  dynaalo  loads  to  anuLll  footlBgs 
has  been  developed*  This  apparatus  Id  designed  to  provide  controlled  loadings 
with  nlnljnua  oo3i?)lexlty.  "Gontrollad"  refers  to  the  ability  to  detonaina  the 
loads  that  are  applied.  Variations  in  the  paranetor  defining  the  loads  are 
introduoedj  but  no  atringont  requirements  are  enforced  on  those  variations 
as  long  as  the  resulting  force-tlae  history  can  be  measured.  The  final  design 
of  tho  apparatus  la  the  result  of  modlfioationa  introduced  during  prellalnary 
tosting  to  achieve  the  desired  end  result.  Thus,  little  could  be  gained  by 
dlseuesing  the  design  development,  and  attontion  will  be  Confined  to  doaoriblng 
tho  present  laboratory  setup.  As  of  this  writing,  the  research  using  the 
dyTiamlc  loading  device  is  preliminary  In  nature  with  the  main  objective  being 
equipment  development,  calibration,  eto.  In  the  eouroo  of  this  work,  however, 
footings  have  been  tested,  loads  measured,  etc'.,  thus  providing  aoaie  quanti¬ 
tative  data. 

Figure  G.l  shows  a  sketch  of  the  oiqporimantal  setup.  Flgurou  0.2 
and  0*3  show  the  dynamic  loading  apparatus  attached  to  the  glass  box  and  sand 

box  respectively.  The  operation  of  the  dynamic  dovlto  is  relatively  slngsle. 

With  the  solenoid-operated  valve  closed,  pressbre  Is  roloased  from  the  nitrogen 
bottle  into  the  air  accumulator  and  hydraulic  accumulator,  so  that  pressure 
builds  up  in  the  hydra  ilic  fluid  behind  the  valve.  To  apply  e  dynamic  load 

the  valve  is  opened,  allowing  hydraulic  pressure  to  be  applied  to  tho  piston. 

The  load  transmitted  to  the  footing  is  measured  by  a  force  waeher^.  The 
equipment  described  herein  produced  rise  times  of  from  2  to  $ma  In  the  loads 
acting  on  the  footings.  The  displacement  of  the  footing  is  msasui’ed  by  a 
linear  vnri'ble  tronsdvio-rr  the  assumption  being  that  this  measurement 
represents  the  displacement  of  the  center  of  the  footing. 

The  e^qjerimentol  studies  with  the  dynamic  apparatus  have  been  pre¬ 
liminary  in  nature.  An  attempt  is  always  made  to  obtain  a  reasonable  amount  of 
data,  but  primary  emphasis  to  date  has  been  on  checking  out  the  equipment. 

The  preliminary  tests  were  conducted,  vdLth  the  displacements  and  loads  being 

^  Lockheed  Electronics  Model  WR7S  High  Sensitivity  Force  Washer, 
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rpcoi'clf'cl  by  a  '''ol'irold  ohota'^rnph  of  thn  oscilloscope  screon.  Althouph  these 
records  were  sufficient  for  this  prpliminnry  rrork,  s  I'onsolidstod  recorder  lias 
been  ^  ntrotiuc*}d  for  experiments  where  the  data  are  of  prino  ^riportanct;. .  Thu 
following  paragraphs  present  a  general  descripti  on  of  the  oparation  of  the 
dynamt c  apparatus, 

Figures  Q.li  and  Q,S  show  two  examples  of  the  records  obtAined  from 
the  oscilloscope.  Figure  Q.h  shows  a  moro-or-less  linearly  decaying  force  with 
a  relatively  small  maxijnum  ^splacement,  .approximately  0.1  In.  Figure  0,5 
s^'owo  a  larger  displacement,  approximately  1,0  in.,  and  a  force  that  drops  off 
after  nbe  t  IHO  ms. 

Figure  0,6  and  0.7  are  examples  of  records  obtained  by  the  Consolidated 
recorder.  These  records  read  from  right  to  left.  Both  displacement  records 
are  beyond  the  linear  range  of  the  recorder  —  this  p:q5lainB  the  apparent  recovery 
after  the  initial  peak.  The  records  aro  meaningful  up  to  the  maximum  displacement, 
but  beyond  that  point  it  is  only  possible  to  deduce  that  the  footings  stopped 
movi  ng,  Fi  mires  0.6  and  0,7  show  respectively  a  relatively  long  and  a  relatively 
short  duration  of  load  application. 

Ihe  advantage  of  the  records  obtained  from  the  j6onaolldated  recorder 
for  quantitative  data  is  obvious.  For  W  s  reason,  it  is  anticipated  that  this 
method  will  be  used  for  the  remainder  of  the  d.vnam’ic  experiments.  By  varying 
the  initial  pressure,  the  status  of  needle  valve,  etc,,  it  is  possible  to 
vary  both  the  shape  and  magnitude  of  the  load  obtained  from  the  system.  Current 
work  is  still  concerned  with  learning  to  control  these  paraiiraters'. 
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Fig.  Q.l  ffiETCH  07  DYNAMIO  APPARATO3 


Fig.  G,2  DPiAHIC  APPARATUS  ON  OUSS  BOX 
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Fig.  G*3  HUJAMIC  APPARATUS  ON  SAUD  BOX 
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(a)  Raoord 


A  fi  M  O  U  B 


(b)  Callbra,t<lon 


Pig.  O.Ii  RECORDS  FOR  Ig-in.x-l-ln.  FOOTING.  TEST  N0.6 


Records  read  right  to  left.  Top  trace  is  force; 
Bottom  trace  is  displaoemenW 
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(b)  Calibration 


Fig.  0,5  RECORDS  FOR  U-in.x-U-in.  FOOTING.  TEST  NO,  2 


Records  read^^  right  to  left.  Top  trace  is  force j 
Uottcm  trace  is  dicpiacement. 
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(a)  Dlcfplftceianot  Roo >rd 


Inegeaainfi  T3j» 
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(b)  Force  Hooord 


IncreaslnR  Time 


Fig.  G.7  .TYPICAL  CONSOLIDATSD  RSCOHDS.  TEST  I'O,  10 
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